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ABSTRACT 

Certain measurements of the electrical and mechanical 
properties of transducers are of importance to designers. 
These can be determined from acourate complex admittance and 
impedance measurements made throughout the frequency range 
of interest. Manual collection of needed data is a 
time-consuming process that is prone no error. The 
computerized system herein described substantially reduces 
the time requirement and produces more accurate output than 
can be obtained utilizing manual, methods. In addition to a 
general discussion of equivalent circuits and the 

instruments employed, this report includes samples of plots 
and calculated parameters for piezoelectric and 
magnetcstr ictive transducers, the experimental comparison 

with traditional manual methods and the program listing for 
two versions of the system allowing varying amounts of 
operator options. 
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I. INTRODUCTION 



Transduction is defined in an extremely general sense as 
any conversion of one form of energy into another. The 
devices responsible for this conversion, in either 
direction, are transducers. For this author, the interest 
is the conversion of electricity into sound and vice versa. 
Commonly, a transducer that converts electrical energy into 
sound is called a projector, loudspeaker, source or 
transmitter. Those serving the reverse function are 
commonly called hydrophones, microphones, or receivers. 
Seme types of transducers are reversible and are capable of 
filling either role. This study is restricted to the 
analysis of reversible transducers. 

A. MECHANISMS OF TRANSDUCTION 

There are two basic types of mechanisms commonly used in 
reversible electroacoustic transducers designed for 
underwater sound, one involving interaction between 
mechanical motions and forces and electrostatic fields and 
the other between mechanical motions and forces and magnetic 
fields. Eoth types are amenable to analysis using 
equivalent electrical circuits; but, due to the basic 
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differences in conventions used to describe forces in 
electric and in magnetic fields, the details are different. 
Magnetically coupled transducers of the moving coil type 
have the diaphram attached to a wire coil suspended in a 
magnetic field. It is the interaction cf the current in the 
coil and the magnetic field that induces a force on the coil 
and causes the diaphram to move and radiate an acoustic 
wave [Ref. 1]. These classifications may be further 
subdivided based on the motor mechanisms involved (i.e. 
electrodynamic, elecor osta tic, magnetic, aagnerostrict ive, 
and piezoelectric) [Ref. 2]. The classification systems for 
transducers vary depending upon the analysis being done. 
Descriptions of various classifications are presented by the 
previous references as well as by Heaslip [Ref. 3], and 
Sherman [Ref. 4 ]. 

B. PURPOSE OF THIS STUDY 

All electroacoustic transducers may be considered as 
being composed of electrical elements so that the analysis 
of the transducer is greatly simplified if an equivalent 
electrical circuit is used. 

While the most common transducers currently in use in 
underwater sound are made of piezoelectric materials, -.he 
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aim of this study is to provide a generalized, 

semiautomat ed means of obtaining meaningful transducer 
measurements on any type of transducer with the operator 
having mininal knowledge of the device. 

The equipment set-up utilizes basic devices (e.g. 
voltmeter, sweep oscillator, Dranetz impedance meter, 

amplifier) interfaced to, and controlled by, a desk-top 

computer. This interfacing provides more precise data 
acquisition than can be obtained from manual plots 

interpreted by an operator. The time saving factor is a 

major motivation for this choice. 

Certain measurements of the electrical and mechanical 
properties of the transducer are of importance to designers. 
In some cases, it is desireable to have the resonant 
frequency fall within a particular range. More likely, 

there will be concern for a high electromechanical coupling 
coefficient. Efficiency nay be the major concern or perhaps 
it is desired to have high losses and, therefore, a low 
mechanical quality factor with consequent "flatness" of 
response [Ref. 5]. 



Direct 


measurement of 


these values 


is not always 


possible. 


However, many of 


the electrical 


component values 
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needed to calculate mechanical parameters can be readily 
obtained from properLy annotated impedance and admittance 
plots for a particular transducer. 

Currently, there are automated systems designed to 
analyze transducers, such as the WQM-12 [Ref. 6]. Few are 
small enough to be readily transported to the transducer as 
is this system. 

C. FORMAT OF THE REPDRT 

The basic physical and electrical theory used in 
preparation of this report is presented in Chapter II. 
While more thorough analyses of transducers and equivalent 
circuits may be found in some of the references, only data 
pertinent to the measurements obtained by this system are 
presented. 

Chapter III expLains the interfacing between various 
components of the equipment set-up and the computer. The 
operator/computer interaction is also addressed. 
Representative samples of computer output are included. 

Chapter IV briefly describes the testing and evaluation 
of the system and program. The results and conclusions are 
presented here. 
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omputer language 



The computer programs written i 
are included in Appendix A. 



n BASIC C 
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II 



THEORY 



Since much of the physics 
discussions of transducers 
electric-coupled and magnetic 
report will make a general 
differences where they occur, 
chapter, contains a listing of 
and refers to symbology used 
cited. 

A. GENERAL 

The approach usually taken 
is to use an electric circuit 
simplest form one may use a tw 
representing the electrical t 
mechanical terminals. Analogs 
voltage and current or vice ve 
transducer. 

The canonical equations de 
two porn transducer may be writ 

V-V 1 + T em*U 

and 



and engineering involved in 
is common to both 
coupled transducers, this 
approach and point out 
Table I, at the end of this 
symbols used in this report 
in some of the references 

in electro acoustic analysis 
o model the transducer. In 
-port network with one port 
rminals and the other the 
for force and speed may be 
sa depending on the type of 

cribing the behavior of the 
e r. as: 

( 2 . 1 ) 
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( 2 . 2 ) 




where V and I are the voltaga and current at the electrical 
port and F and U are the force and velocity at the 
mechanical port. I em is a transaction coefficient relating 
the electromotive force at the electrical port to the 
velocity in the mechanical net, and T me is a transduction 
coefficient that relates the force developed at the 
mechanical port of a two-port network to the current in the 
electrical mesh [Ref. 8]. <t> is the transformation factor 

and may be defined as 4>=T ejn /Z 0 . The canonical equations for 
a magnetically coupled transducer of the moving coil type 
are: 



where 3 is the magnetic field and 1 is the length of the 



electrodynamic transducers (moving coil) as written using 



F= - B • 1- I + Z m *D 



(2.3) 



7=Z 0 - I + B-1-0 



(2.U) 



wire in the coil 



(Note: This choice is valid for 



B •l=T em . Magnetostri ctive transducers obey the same set of 



equations but T em £ B-i.) 
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For electromechanical coupling. Figure 2.1 shows typical 



circuits and transformations ^Ref. 7], 




. l 9m 



Vr^HSIH 0 



in = 



Fig. 2.1. Typical Circuits and Transformations 



In the vicinity of resonance, which is of concern here, 
rhe Mason circuit or modified Van Dyke circuit is valid. It 
is also valid for low frequency values which will also be 



needed 


for some 


calculations * Ref. 8]. 


This circui 


shown 


in Figure 


2.2 where L m , R m , 


and C m repr 


mechan 


ical values. 


: Q is the clamped or 


blocked or 



capacitance . 
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o 




Figure 2.2. (a). Modified Van Dyke or Mason circuit (b) . 

Basic Transducer Circuit. 

Kinsler, and others [Ref. 1 ] r define the following 
measurable mechanical and electrical impedances of a system: 

Blocked electrical impedance (ohms) =Z Q = V/I I u = q _ 

Free electrical impedance (ohms) =Z =V/I|_ 

F = 0 (2.64 



Open 



circuit mechanical impedance 



(N s/m) =Z m =F/U( I=0 

(2.7) 
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Short circuit mechanical impedance (S • s/m) = 2 mm =F /U I v _q 

( 2 . 8 ) 

Magnetically coupled transducers may be represented by 
the circuit of Figure 2.3, where ^ is a transformation 
factor. 4> m = B • 1 for moving coil transducers [Bef. 2], 




o 







o 



Figure 2.3. Equivalent Circuit for a Magnetically Coupled 
Transduce r. 

Important properties of the transducer may be found by 
studying the electrical impedance of the system as a 
function of driving frequency. The unloaded driving point 
electrical impedance may be determined from the canonical 
equations (2. 1-2.4) and is given by 

Z * ?/I| F.(T Z o ‘ '•Wme'Zml <2-51 
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Similarly, the short circuit mechanical driving point 

impedance is given by 

^mnT^^lv=0 (2.10) 

The electrical impedance is influenced by the motion of the 
coupled mechanical system, as Indicated by the second rerm 
on the right hand side. This term is referred to as the 
motional impedance, and defined by Hunt [Bef. 2], as 

^mot = (“^em ‘^me^ (2.11) 

Similarly, the mechanical driving point impedance is 

influenced by the electrical load. 

B. USE OF COMPLEX IMPEDANCE /A OMITTANCE PLOTS 
1 . The Resonance Circle 

Graphical displays of impedance and admittance data 
can be extremely useful in the analysis of a transducer. 
When properly scaled and annotated, quantitative results may 
be obtained from the diagrams. Generally, in the vicinity 
of resonance, the eLement is regarded as having a single 
degree of freedom and the locus of points of the complex 
impedance is a circle in the neighborhood of the resonance 
frequency [Ref. 7], [Ref. 3]. Figure 2.4 and Figure 2.5 
are typical resonance circle diagrams where, in Figure 2.4, 
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the electrical reactance is plotted as a function of 
electrical resistance with frequency as a 
parameter (Bef. 8]. F4 and f5 are the frequencies of the 
half power points and f8 is the resonance frequency. Dashed 
lines show the same transducer under load conditions. 





Fig. 2.4. Complex Impedance Diagram 



(Note: Reactance may be entirely negative for an 

electrically coupled transducer! . 
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Fig. 2.5. Complex Admittance Diagram 



In Figure 2.5, f2 and f3 equate to the frequencies 
of the half power points. FO is the frequency of resonance. 
The dashed line would be the same ideal transducer tested 
under load conditions. (Note: Susceptance may be entirely 
negative for a magnetically coupled transducer.) 

Motional impedance components may be obtained using 
the following formulae: 
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R mot ” R ° (2 . 12) 

x mot =X E + ^L 0 (2.13) 

Wx*h Z mo j. - R mo t + 3 ^mot (2.14) 

i 

The motional admittance components are given by: 

G mot * g E - (’/Bo) (2.15) 



^mot = Sr, - j C t 



(2 . 16) 



with Y mo t = J mot ~ j ®mot 



(2.17) 



The frequency of maximum power output at constant voltage 
input is called the mechanical resonance end is represented 
in Figure 2.6 by f8 [Ref. 9 ]• 

In the motional admittance diagram. Figure 2.7, fO 
is the frequency of maximum conductance and resonance. F2 
and f3 are the half power points. 

2. Component Responses 

Impedance measurements are made at the electrical 
terminals of the transducer. The behavior of the real and 
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Fig. 2.6. Motional 

imaginary components of the 
observed as the frequency 
sufficiently slow rate so 
quasi-steady-state. When the 
so that it cannot move, the bl 
be measured and used to calcul 
it is difficult to adequately 
frequency spectrum, values m 
below resonance and tie interm 



Impedance Diagram 



complex 


impedance may 


be 


spectrum 


is 


swept at 


a 


that 


4- U - 

w U ^ 


system 


is 



mechanical system is clamped 
eked electrical impedance may 
te motional impedance. Since 
lamp a device over the entire 
y be obtained far above and 
diate values inferred. This 
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Fig. 2.7. Motional admittance Diagram 



is especially important for 
curves whose clamped value 
frequency than do the real 
impedance. 

The maximum value of 
to the mechanical resistanc 
circle. This value will be o 
efficiency of the system. 



the reactance and susceptance 
vary more drastically with 
components a f admittance or 

actional impedance will equate 
and is the diameter of the 
importance in calculating the 
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in put 



Plots of various input electrical parameters 
(conductance (G E ) and susceptaace (B E ) or resistance (R E ) 
and reactance (X E ) versus frequency) can provide all 
of the information necessary to obtain the electrical and 
mechanical quality factors for the system. These values 
indicate the sharpness of the resonance, 
a. Impedance 

The circuit of Figure 2.3 is useful for analysis 
of the behavior in the vicinity of resonance. 




Fig. 2.8. Equivalent Circuit Diagram 



A plot of the input electrical impedance components versus 
frequency will resemble Figure 2.9. 
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o xi 0 W 



Fig. 




2.9. Resistance and Reactance Plotted Versus Frequency 



b. Admittance 



The motional admittance equations are 



Y mot = J mot + j 3 mot 



J ^mot~j X mot_ 



R mot + j X mot H mot"“j Xm ot 



®mot 



^mot 



Vot + ^mot 



®mot + ‘^mot 



( 2 . 18 ) 
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Plots of conductance and susceptance versus frequency for 
Figure 2.8 will resemble Figure 2.10. The frequencies where 
the susceptance (B) is zero correspond to the electrical 
resonance and antiresonance frequencies. 




Figure 2.10. Conductance and Susceptance Plotted Versus 
Frequeuc y 

C. EVALUATION OF TRANSDUCER PARAMETERS 

This system is designed to cbtain useful information on 
devices used in transmitting and receiving in water. As 
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such, it is first necessary to obtain data in an unloaded 
situation in order to be able to separate the effects of 
acoustic loading. For transducers designed for use in 
water, air normally provides an adequate medium for 
evaluation cf the acoustically unloaded properties. Placing 
the transducer in water means the two-port network will now 
reflect the sum of the mechanical and the load impedances. 
Then, according to Hunt [fief. 2], 

Z T = Z m + Z 1 = X,) (2.19) 

or 



z ra + Z 1 = < R m f R l) P + j ’ 2 • Q • p) 



( 2 . 20 ) 



where p is a frequency par ameter-- (p=. 5 ( (u4j 0 ) - (u 0 /u) ) and Q 
is the quality factor measured for the resonance (i.e. 
electrical or mechanical depending on the type cf 
measurements) . 

Data obtainable from admittance measurements will 
provide needed valnes for determining the motional 
admittance circle. The maximum value of conductance will 
occur at the frequency of elactical resonance. The qua lity 
factor can be determined from the frequencies an which the 
conductance assumes half the diffferer.ca of its maximum 
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value and its blocked value (i.e. Q = f /(f u -fj). These 
frequencies will also coincide with the frequencies (f2 and 
f 3 ) of the local maximum and minimum susceptance values near 
resonance. Blocked capacitance may be calculated from 
measurements of susceptance far above and/or below 
resonance. For an electricaLly coupled transducer with 
small dielectric losses, valid measurements of C 0 may be 
made very far below resonance. For magnetic coupling, 
measurements must be taken far above and far below resonance 
and the values between interpolated linearly between these 
two. 

Impedance data will provide the frequency of 
antiresonance. The frequencies of the half values for the 
difference between the maximum resistance and the blocked 
resistance will correspond to tae maximum and minimum values 
of reactance and allow for calculation of a mechanical 
quality factor. Blocked capacitance or inductance 
(depending on the type of transducer) may be calculated from 
reactance data obtained far above and/or below resonance. 

The following paragraphs briefly describe the definition 
and calculation of various constants. 
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1 



• Qua lity Fa ctor = 2 

This factor indicates the sharpness of resonance. 



Q = 



frequency of max. G 

freq. difference of 
and max. B or X 



or R 
min. 



f (resonance) 
f (upper) - f (lower) 



( 2 . 21 ) 



2. Electromechanical coupling coefficient = (ef f ) 

This coefficient indicates how tightly the 

electrical and mechanical meshes are coupled [Ref. 2]. 

K 2 (eff) = 1-(f3/f0) 2 for electromagnetic coupling 

( 2 . 22 ) 

and 

K 2 (ef f) = 1-(fD/f8) 2 for electrostatic coupling 

(2 .23) 

where f8 is used to indicate the frequency of maximum 
resistance and fD indicates the frequency of maximum 
conductance . 

3. Static Coupling. Coefficient = K£ 

This term is the ratio of the stored mechanical 
energy to the total energy stored in the device [Ref. 1 ], 
[Ref. 5]. 
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(2. 24) 




Co + C i 



4. Shunt C apacft a nee = C Q 

This is the blocked capacitance of the system. 
There is no simple method of measuring this quantity (with 
the possible exception of Instruments designed for a 
specific frequency). !i easurements must be made at a 



frequency 


far 


re mov ed 


from frequencies 


of mechanical 


resonance 


in 


order 


to minimize the 


effects of 



electromechanical coupling. If the capacitance ratio 
(C^/C 0 ) is greater than fifty, the inaccuracy should be 
minimal and may be considered negligible. 

1 3 

C 0 where <*> << co Q (2.25) 

coX q co 

5. Cl a raped Resistance and Reactance = Rq ana Xq 

R q and X 0 may be closely approximated by obtaining 
data far below and above resonance and interpolating between 
since there is a positive frequency dependence. However, 
for many practical purposes, (and for electrically coupled 
transducers), data ma/ be obtained far below resonance which 
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provides a reasonable approximation to the actual value 
throughout the frequency range of interest. These values 
are assumed to represent those which would be measured if 
the mechanical parts of the system were blocked or clamped 
to prevent motion. 

6. Hocked Inductance = Lq 

This is the inductance of the blocked system 
(usually with magnetic coupling) and is obtainable from data 
collected far below and above resonance. 

X(interpola ted) 

L = (2.26) 

“o 

7. Capacitance Ratio = 12! 

This is the ratio of the blocked capacitance (C Q ) to 
the motional capacitance (C,). 

C Q PI ( 1 — K 2 ) wher<= K* is the 

'R' = static coupling 

C ! 8 :< 2 coefficient 

(2.27) 
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8. Fi g ur e of Merit 



The 


figure of merit 


i s 


little 


mentioned 


in the 


references 


probably because 


it 


is an 


arbitrary 


for mula 



depending upon a designer's criteria. Miller [Ref. 5], 
mentions K 2 -Q m as one possibility and utilization of the 
coupling co eff icient as another. No specific calculation of 
this value is done ia this report, leaving it to the user to 
insert as desired. 

9. Pi a meters 

Impedance or admittance circles ace transformed to 
motional diagrams by subtracting blocked values or blocked 
values multiplied by angular frequency as discussed in 
section B. 1 above. The diameters of the resulting circles 
will provide needed values for efficiency calculations. 
and signify the diameters of the motional circle for the 
transducer in water and in air, respectively. 

D. EFFICI2 NCI 

The overall efficiency of a transducer is the ratio of 
the power delivered to an external load connected at the 
output terminals to the total power at the input. 

At resonance, this can be easily calculated from known 
values. For an eLectrically coupled transducer, the 
efficiency is given by: 
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E 



(2 .28) 



G D~ 



where and are the diameters of the motional admittance 
circles and where G is the conductance at resonance. 

Similarly, the efficiency at resonance for a 
magnetically coupled transducer is given by: 

E = (2.29) 

R-D* 

where these are diameters of the motional impedance circles 
and R is the resistance at resonance [Ref. 10]. 

The mechanical power utilization factor is part of the 
overall efficiency. It is given by R l / which can be 

expressed in terms of loaded and unleaded values of Q or by 
the diameters of circLes. 



' m 




(2.30) 



(Note: All values should be taken from either admittance 

data or impedance data but should not be mixed.) 
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Once the efficiency at resonance is determined, 
utilization of the impedance data of a magnetically coupled 
transducer allows for calculation of the lagging phase angle 
(( 3 ) • 



P = 0.5* ArcCos ((R res -R 0 > /D w ) (2.31) 

The frequency parameter for maximum efficiency may be 
calculated as 

p= {D w Sin(2-p)| /(4-R q . Q w i (2.32) 

and from the earlier relationship between d and the 
frequency ratios, the frequency for maximum efficiency f or a 
magnetically coupled transducer may be calculated. 

f = fO ( 2 • p + <p* + 1) vfe) (2.33) 

The frequency of maximum efficiency may not coincide 
with either the resonant frequency or the frequency of 
maximum admittance. For the magnetostrictive transducer, 
this frequency of maximum efficiency is usually greater than 
that at resonance. 
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A development oy Camp * Ref . 8], shows that a 

magnetostr ictive transducer may produce more acoustic power 
at the resonant frequency even though efficiency may be 
less. Normally, the frequency for optimum acoustic power 
will fall between the mechanical resonance frequency and the 
frequency of optimum efficiency. This is due to eddy 
current losses and magnetic hysteresis. 

Since transduction losses are low, maximum power output 
for piezoelectric devices is limited mainly by voltage 
limits. The optimum operating frequency is accepted as the 
frequency of maximum admittance * Ref . 8]. 

If one assumes that the transducer has a high Q (sharp 
resonance) and that the blocked resistance is constant 
throughout the frequency range for the impedance circle, 
then the potential efficiency may be expressed as 

Pot. Eff. = n --— (2.34) 

(Hmax) 1 ^ + (5 min ) 1/2 

The potential efficiency represents the maximum efficiency 
available for the most favorable loading conditions. This 
equation deletes the requirement for data from a loaded 
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condition but requires the air motional impedance loop to 
very closely approximate a circle [Ref. 9]. "High Q" as 
used here implies values greater than thirty. If the air 
motional impedance loop is not approximately circular and/or 
the Q is low, then, from reference 2, 



Pot. Eff. 



(1 COS 2 P ( V* - (1-j& SIN 2 (3 ,V» 

(1 «■ COS 2 (3 l V* + (1- SIN 2 p ) ^ 

(2 .35) 



E. SUMMARY 

There is helpful information to be obtained from both 
impedance and admittance data. Some transducers more 
readily reveal properties through one type of diagram than 
another. Normally, admittance data is collected for 
electrically coupled systems, while impedance data is better 
for studying magnetically coupled systems [Ref. 2]. Due to 
the ease of obtaining data with the system proposed hers, 
both sets of data and diagrams nay be rapidly obtained. 
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TABLE I CCONT) 
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III. ADAPTATION TO A ZOMPDTERIZSD SYSTEM 

In-phase and out-of-phase components of the admittance 
or impedance are needed in the analysis of a transducer. An 
initial attempt was made to use a lock-in analyzer, but 
ultimately a Dranetz impedance meter was used. 

A. EQUIPMENT SET-OP 

Data needed foe the calculation of properties of 

transducers may be obtained using a lock-in analyzer. The 
magnitude and phase or in-phase and quadrature components 
will indicate resonances aid provide the necessary 
information . 

The initial set-up utilized the Princeton Applied 

Research Company Model 5204 Lock-in Analyzer. Phase 

information would have allowed for desired calculations. In 
experimentation, the circuit required for obtaining valid 
phase information became too complex to be practical. The 
lock-in analyzer is designed to operate in the 1 Hz to 100 
KHz frequency range although it was determined to be 

accurate at higher frequencies. At low frequencies, there 
is a low source resistance and the signal to noise ratio 
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became too small to provide a reference. It was concluded 
that complex impedance measure n ants would be another means 
of getting the necessary data. 

As briefly mentioned in the introduction, there are many 
methods for obtaining impedance information. Since this 
au+hor desired a portable system usable for underwater 
transducers, a bridge network seemed reasonable. 

The set-up used to design tnis system is shown in Figure 
3.1. The Hewlett-Packard HP-35 computer directs the 

peripherals for the data taking. The printer and plotter 
are non-essential to the system and were used to provide 
larger format graphics and printout: than are available with 

the built in thermal printer provided with the HP-85. 

1. Hewlett-Packard HP -85 Computer 

The HP-85 is an eight bit microcomputer that 

utilizes BASIC computer language. It has a 127 millimeter 
diagonal black and white electromagnetic-deflection CRT. A 
32 character per line thermal printer/plotter is part of the 
unit. The computer has 16K bytes of read/write memory of 
which 14,579 are available to the operator. For this 
program a memory module is needed to expand the memory to 
32K bytes. Programs or data may be stored on, or read from. 
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magnetic tape cartridges. To interface with peripheral 
equipment, an I/O ROM and an interface card were added to 
provide HP-IB (IE3E-438) instrumentation interface 
capabilities. 

2. Synthesizer /Function Generator 

The Hewlett-Packard Model 3325A Synthesizer/Function 
Generator can produce five different continuous waveforms. 
For this application a sine wave was desired for which this 
instrument has a frequency range of 1 microhertz to 20 
megahertz. Frequency may be specified with up to eleven 
digits of resolution. Output amplitude is 1 millivolt to 10 
volts paak-to-peak into a fifty ohm load. This model is 
fully programmable through the rear panel Hewlett-Packard 
Interface Bus (HP-IBl . 
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Figure 3. 1 Equipment Set-up 



current is set by the range switch. The voltage developed 
across the unknown impedance is amplified and fed to two 
resolvers. For measuring admittance, a constant voltage is 
applied to the unknown. The current through the unknown is 
measured by means of voltage developed across a small 
range- selected resistor placed in series with the unknown. 
The resolver outputs are proportional to conductance and 
susceptance of the n known. The claimed accuracy of the 
instrument is plus or minus two percent full scale 
amplitude. The Dranetz meter requires an input signal 

voltage of 3 Vrms. In taking measurements on an unknown 
transducer, initially the frequency of interest may not be 
known as well as the appropriate scale settings for the 
meter. Measurements made on the Dranetz meter using a 
variable resistance box and voltmeters indicate linearity 
throughout the range settings. However, it is necessary to 
insure measurements do not cause either mater range to be 
exceeded. The limit indicator lights are activated when the 
unknown impedance or admittance is substantially greater 
than the full scale setting. fhen this occurs, an internal 



relay shunts the meters and the DC output terminals 



to zero. 
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• Da ta Acquis it i. on/ControL (Jn it 

The Hewlett-Packard Model 3497A Data 
Acquis it ion /Control Onit is used here as a voltmeter and 
scanner to measure the DC output voltage proportional to the 
deflection of the meters of the Dranetz. The instrument 
measures and displays voltage values to five and one naif 
digits. It is internally triggered by software command 
during the data taking but takes the next reading as soon as 
the first is completed for obtaining averages. 

5. Pow er Amplifier 

This instrument is not essential to the operation of 
the system. The Dranetz meter is designed to accept direct 
ac voltage input; however, the particular meter used 
exhibited more stability when the amplifier was included in 
the set-up. The amplifier had unity gain and was only used 
to lower the effective output impedance of the frequency 
synthesizer from fifty ohms to two milliohms. 

6. General Discussion 

This equipment set-up is suitable for an initial 
system. A system designed for regular and rcurine use 
should include a different impedance measuring system than 
the Dranetz. Because of the design feature of the Dranetz, 
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a great deal of operator interface is required to change the 
various different scale and raige factors which cannot be 
set directly by the computer. It is too easy for an 
operator to miss changing one setting or to overlap 

frequency range scales in the spectrum sweep and, thereby, 
obtain erroneous data. Also, the scale factors must be 
manually entered prior to the computer doing any cal- 
culations. The Dranetz also had annoying features, such as 
drift of the zero, excessive noise on low impedance ranges, 
and DC output at full scale deflection of 0.8 - 0.95 volts. 

B. DISCUSSION OF THE PROS RAM 

The HP-85 is the controlling unit directing the other 
instruments to take or send data as needed. Memory storage 
is handled by this microcomputer. 

The goal was to design a program than would minimize the 
time necessary to collect pertinent data while maximizing 
accuracy. The program was to accommodate different media, 
allow for collection of admittance or impedance data or 
both, and work for transducers with either electric or 
magnetic coupling, while being easy to operate. This was 
accomplished; but due to the large differences in properties 
of different transducers, a second program was written tc be 
used by more experienced operators. 
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One program is designed so that a person with minimal 
knowledge of the transducer or measurement procedures may 
run it and obtain meaningful information. It is suggested 
that the operator be familiar with the operation of the 
interfacing instruments (especially the Dranetz meter). An 
equivalent program is designed so that the operator may 
select the bandwidth for data collection and may specify the 
frequency for measurement of tie blocked values (done as a 
percentage of the resonant frequency) . 

The program is broken into three sections. The first 
gives the operator an overview of the conductance and 
resistance of the transducer uader test. Part two is the 
data collection for desired measurements with options for 
plots and lists of the data. The final section performs the 
calculations for motional data, efficiencies, and other 
desired outputs. The folloviag subdivisions address the 
program in detail. Figure 3.2 is a flow chart of the 
program. The programs written in BASIC are included as 
Appendix A. 

1. The Search 

This first portion has been included to allow for 
visual and graphic display of tie response of ohe transducer 
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over a selected frequency range. The operator is required 
to insert information on the type of coupling (magnetic or 
electric) of the transducer and the medium (air or water) in 



which measurements 


at re 


to be taken. 


The operator is 


instructed to 


set 


t he 


Dranetz for the 


measurement of 


admittance data 


and 


t o 


specify the range 


to be covered in 



the frequency sweep. The operator also enters the desired 
voltage to be sent to the Dranetz. This manual voltage 
input is necessary since at higher frequencies lower rms 
voltages are necessary to keep the ac voltage meter on the 
Dranetz at 3.0 7. 

Through the Data Acquisition Unit, the relative 
amplitude of the real part of the admittance is sampled at 
300 equally spaced points throughout the selected frequency 
range. A graph is displayed on the CRT and then copied on 
the thermal printer or the HP-33 for later reference. In the 
event that some parameter needs to be changed, the operator 
is given the opportunity to rerun this portion. 

The operator is instructed to specify a relative 
amplitude (from 0 to 1.2) that is less than the maximum of 
the admittance peak of interest. The computer prints the 
relative a nplitudes greater than this value and the 
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associated frequencies. This allows the operator to select 
the resonance of interest and provide the center frequency 
and halfwidth to be considered in finding the resonant 
frequency and quality factor (2> . A comparison is done to 
find the absolute maximum value for amplitude and the half 
power points for the determination of 2. The resultant 
factors are printed and the option to rerun is again 
extended. 

This same subroutine is again used to obtain similar 
data for the real part of the impedance. The subroutine is 
an adaptation of a part of a procedure designed by 
Conte [ Hef . 13 ]. 

Although two repeats of the subroutine take about 10 
minutes, the operator is given sufficient information on one 
impedance and the admittance responses of the transducer to 
decide which type of measurements will provide the most 
accurate information. 

2. Data Collection 

The computer has the frequency synthesizer sand a 
frequency far below resonance. The meters of the Dranetz 
are zeroed by the operator and set for full scale 
deflection. An average of ter. measurements for each 
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Figure 3. 2 Flow Chart of the Program 
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Figura 3. 2 (Continued) 
56 



terminal is placed la memory to be used as normalization 
factors for ail later DC output voltage measurements. 

While still at this frequency far below resonance, 
an average of ten measurements will be made, normalized and 
scaled (using operator input for the scale factors from the 
Dranetz instrument) and stored as blocked values of 
conductance, susceptan ce , resistance and reactance. These 
will be used later in computations after being assigned to 
different variables depending on the medium. This one 
measurement is all that is necessary for the electrically 
coupled transducers with small dielectric losses. For the 
magnetically coupled transducers, the values are taken far 
below resonance and stored. The program loops back through 
this section after sending a message to the frequency 
synthesizer to go an equal frequency distance above 
resonance. An average of the two values is stored for later 
use. This averaging is necessary since the eddy current and 
hysteresis losses cause a steady increase in the blocked 
impedance values. 

The decision of tie operator on the type of data to 
be collected is entered. The scale factor from the Dranetz 
instrument is entered and the oomputer goes to a subroutine 
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The 



for data collection. The computer draws the appropriate 
center frequency from the data collected during the search 
and bases the bandwidth for the sweep on the quality factor 
measured previously. Fifty frequencies and the associated 
real and imaginary data are measured and stored from ten 
times the half power bandwidth below center frequency to two 
times the bandwidth below resonance. Two hundred 
measurements are taken in the vicinity of resonance and 
another fifty in a region from two to ten times the 
bandwidth above resonance. (Note: This program reduces the 
range of interest to plus and minus five times the bandwidth 
around resonance when tests are done under load (water) 
conditions. The second program takes data at three hundred 
equally spaced frequencies within the bandwidth specified by 
the operator.) Once the data nave been collected, another 
subroutine is entered to find the minimum and maximum values 
for the two sets of data and their associated frequencies. 
This is accomplished through a comparison to a minimum (or 
maximum) value until a smaller (larger) one is found. A 
reassignment is done and the comparison continues throughout 
the data. The maximum value and che two adjacent values are 
then fitted to a quadratic equation which is used to produce 
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a better estimate of the maximum ani its associated 



frequency. In the event of an error in data collection, 
there is an opportunity to retalce the data before proceeding 
in the program. 

A subroutine provides for a listing of the collected 
data (100 points) in the vicinity of resonance. All of the 
data could be printed; however, the data near resonance 
provides the important information and printing is an option 
open to the operator. 

Plots of the data may be made. The operator may 
chose to make plots of real or imaginary amplitudes versus 
freguency or the plots of real versus imaginary data over 
frequency, or all, or none. Each plot is drawn using a 
subroutine and labeled with pertinent scaling and captions. 
At the end of this chapter. Figures 3.3 through 3.14 are 
examples of plots for a piezoelectric transducer. (This 

example is for a Type 100-5353-010, Serial No. 1213, 

designed for use in an active sonobuoy.) Figures 3.15 
through 3.22 are examples of plots for a magnetostrictive 
transducer. (This example is a Type CtlT-10255, Model QGB, 
Serial No. 318.) 
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3. Da t a Manipula t ion 

In the calculation of notional data, the type of 
coupling of the transducer becones important, 
a. Magnetic Coupling 

For a magnetically coupled transducer, the 
motional impedance data is calculated by subtracting the 
blocked resistance from each resistance measurement. 
Blocked inductance at resonance is calculated using the 
averaged blocked reactance value. The angular frequency 
multiplied by the blocked inductance is then subtracted from 
each value of reactance. In the event admittance data has 
been collected, the motional conductance is obtained by 
subtracting the blocked conductance, while motional 
susceptance may be obtained by subtracting the blocked 
capacitance multiplied by the angular frequency from each 
value. A list of motional data may be obtained. A plot may 
be elected and drawn and labeled with axis scaling done 
using minimum and maximum values of the calculated data. 

Computations are done using the collected data. 
The dynamic electee mechanical coupling coefficient is 
calculated from the resonant frequencies. A new quality 

factor is calculated based cn the frequencies of maximum and 
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minimum amplitudes of the imaginary part of the data 
collected. Values of some blocked measurements, quality 
factors, the coupling coefficient, and resonances are output 
tc the printer. 

b. Electric Coupling 

For an electrically coupled transducer, 
dielectric losses are small and it is acceptible to use the 
values measured far below resonance as the blocked 
measurements. Motional resistance is obtained by 

subtracting the blocked resistance from each measurement. 
The addition of the angular frequency multiplied by the 
calculated blocked inductance to the reactance measurements 
will provide the motional reactance. The data may be listed 
and a plot made. An example of the motional impedance plot 
can be seen in Figure 3.9. 

Motional admittance data may be obtained by 
subtracting the reciprocal of the blocked resistance from 
each conductance value. The angular frequency multiplied by 
the blocked capacitance is subtracted from each susceptance 
measurement. Figure 3.10 shows a motional admittance plot 
for a piezoelectric transducer in air. Calculations similar 
to those for magnetic coupling are done and desired values 
printed. 



5 1 



After measurements are taken with the transducer 



under load (in water) calculations are done to find the 
efficiency and optimal operating frequency (for magnetic 
coupling). This information is output to the printer. An 
opportunity to repeat data collection in either medium is 
afforded before data on air and water temperatures and 
transducer identification numbers are requested prior to the 
end of the program. A sample output (for the 

magretostr ictive transducer Type CMT-10255, Serial No. 313, 
Model QGB) is included as Figure 3.23. 
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hiqure 3.3 Conductance Spectrum for a Piezoelectric Trancducer 
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Fiqure 3.4 Susceptance Spectrum for a Piezoelectric Transducer 
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Fiqure 3.5 Input Electrical Admittance Plot for a Piezoelectric Transducer 
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Fiqura 3.7 Reactance Spectrum for a Piezoelectric Transducer 
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Figure 3.8 Electrical Input Impedance Plot for a Piezoelectric Transducer 
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Fiqure 3.9 Motional Impedance Plot for a Piezoelectric Transducer 
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Figure 3. 10 Motional Admittance Plot for a Piezoelectric Transducer 
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Fiqure 3.11 Conductance Spectrum for a Piezoelectric Transducer 
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Fiqure 3. 13 Input Electrical Admittance Plot for a Piezoelectric Transducer 
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FREQUENCY <HZ> 

-54 1 R VS. FREQUENCY IN AIR 

Fiqure 3. 15 Resistance Spectrum for a Maqnetostrictive Transducer 
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FREQUENCY (HD 
X VS. FREQUENCY IN MR 

Fiqure 3.16 Reactance Spectrum for a Maqnetostrictive Transducer 
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-38 RESISTANCE (0U4S) JL 

MOTIONAL IMPEDANCE PLOT IN AIR 

Fiqure 3. 18 Motional Impedance Plot for a Maqnetostrictive Transducer 
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Flqure 3. 19 Resistance Spectrum for a Maqnetostrictive Transducer 
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X VS. FREQUENCY IN WATER 

Figure 3.20 Reactance Spectrum for a Maqnetostricti ve Transducer 
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Figure 3.22 Motional Impedance Plot for a Magnetos trictive Transducer 



VALUES MEASURED IN AIR FROM ADMITTANCE DATA 

BLOCKED INDUCTANCE = 8.86E-4 HENRIES AT 23956 HZ 
ELECTRICAL QUALITY FACTOR = 189.73 
MECHANICAL QUALITY FACTOR = 29.426 
ELECTRICAL RESONANCE = 23956.0379323 HZ 
MECHANICAL RESONANCE = 23920 HZ 

DYNAMIC ELECTROMECHANICAL COUPLING COEFFICIENT = .003 
BLOCKED RESISTANCE = 33.724 OHMS 



VALUES FOR MEASUREMENTS IN WATER FROM ADMITTANCE DATA 

BLOCKED INDUCTANCE = 9.01E-4 HENRIES AT 23932.48 HZ 
ELECTRICAL QUALITY FACTOR - 18.513 
MECHANICAL QUALITY FACTOR = 19.688 
ELECTRICAL RESONANCE = 23932.4882297 HZ 
MECHANICAL RESONANCE = 23892 HZ 

DYNAMIC ELECTROMECHANICAL COUPLING COEFFICIENT = .00338 
POTENTIAL EFFICIENCY = .20398 
FREQUENCY OF OPTIMUM EFFICIENCY » 24246.61 HZ 
MECHANICAL POWER UTILIZATION FACTOR =* .14645 



Figure 3.23 Sample Computer Output 
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IV 



TEST AND EVALUATION 



Although many of the calculations mala by the computer 
are basic, it was desireable to have a set of data collected 
manually to ensure accuracy of the semiautomated system. 

A. TESTING 

1 . Met hod 

Data were collected for four different transducers, 
all of which were later used to test the computerized 
system. A similar equipment set-up was used with the 
computer replaced by an operator and an X-Y plotter, and a 
Wavetek pulse/functl on generator in place of the 

synthesizer/function generator. 

The frequency spectrum was swept manually with plots 
made of the real and imaginary components versus frequency. 
From these plots, amplitudes and frequency values were 
obtained and calculations similar to those done by the 
computer performed. 

2. Com pari s on s 

For the magnetostrictive transducer used here, the 
data output of the computer corresponds very closely to that 
collected manually. The blocked inductance was found to be 
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the same calculated by either method to five digits. 
Resonance frequency variations up to. ten hertz may be 
attributed to temperature variations; and these frequency 
variations will affect the coupling coefficient. The 
efficiencies compared to within five percent. 

Comparisons of the piezoelectric data ware as 
favorable as that for the magaet ostrict ive transducer. The 
blocked capacitances varied by up to six percent ana there 
were small resonance frequency differences in some cases. 

B. EVALUATION 

The resonance frequency variations noted may be 
attributed aainly to temperature changes from one experiment 
to another. Larger deviations noted in a few cases are due 
to inaccuracies in reading graphs or in determining a true 
maximum for a very low 3 device as is the case for many 
transducers when placed under load. For the examples used 
here, the computer data are considered more accurate since 
it measures and compares data to five digits, far better 
than the average eye can do with unsealed data. In the 
computerized system, the data collected in the vicinity of 
resonance are compared to find the absolute maximum value 
and the associated frequency. This point and two adjacent 
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points ars then used for a parabolic fit to find the precise 
frequency for maximum component amplitude. This provides 
accuracy in the resonance frequency far in excess of what 
may be achieved manually by an operator. 

A 'delay time' to allow stablization of the system 
(dominated by the response time of the Dranetz meter) was 
determined to be one hundred milliseconds. In manually 
sweeping the frequency spectrum, stablization time is 
difficult to achieve. This may result in the appearance of 
spikes in the plots which increase the difficulty of manual 
evaluation of the collected data. Once this stablization 
time was determined, all wait times in the computer program 
were ensured to be greater than this value with the longest 
waiting period during determination of the shunt and 
normalization values. The second longest wait period occurs 
during the data collection. Here, a period of two hundred 
fifty milliseconds has been used in the vicinity of 
resonance to minimize the chance of error. The added minute 
for data collection that this additional wait time requires 
is of minor inconvenience compared to the improvement in 
accuracy. 
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Transducers may exhibit a moderately high Q in air 
(above 25) . For these (and most resonances with a Q above 
5) the program with set frequency ranges for data taking 
works well. However, these same transducers may have very 
low quality factors under load (in water). To avoid data 
review over a frequency range beyond the upper limits of the 
instruments and perhaps overlapping other resonances of the 
transducer, it is recommended that the program allowing for 
operator input of the frequency range be utilized. 

C. CO NC LOS IONS 

It was noted that most transducers have more than one 
resonance. To be a useful analytical tool, the operator of 
this system needs to know the primary design frequency of 
operation for the test transducer. Should the resonance at 
this frequency cause the admittance or impedance to exceed 
the limits cf the Dra netz meter, the operator may opt for 
data taking of the type to fall within limits. 

Manual collection of impedance and admittance data for a 
single transducer in two meiia (air and water) and 
calculations similar to those done by the computerized 
system utilizing data at three hundred frequencies requires 
on the order of fifteen hours. This system carries out 
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these operations in seventy-four minutes. For most 
transducers, only one type of data (admittance or impedance) 
is normally collected. A typical measurement which produces 
results and the plots similar to Figures 3.15 through 3.23 
would take fifty-two minutes. (Since each collection of 
data allows for four plots and two lists of data, the time 
to run the program collecting impedance and admittance data 
in both media could be reduced by forty minutes and 
collection of one type of data in both media by twenty 
minutes if no plots or lists were made and only the listing 
of transducer parameters similar to Figure 3,23 desired.) 

The major time saving factor for this system would be 
found through the replacement of the Dranetz Impedance M eter 
with an instrument that does not require as much operator 
interface. However, semi-automation as provided in this 
system constitutes a major time-saving over manual means 
even when the equipment set-up is not altered in a major 
way. 

The plots produced using an X-Y plotter and sweep 
oscillator enable measurements of frequency with an accuracy 
of about one percent. Using the computer, the precision is 
better than one ia ten taousar.d. Similarly, for 
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calculations using real and imaginary components of zhe 
impedance or admittance, computer precision is better than 
one in ten thousand as opposed to about one in ten for 
manual manipulations. However, the measurement accuracy is 
at best plus or minus two percent for the Dranetz meter. 
The advantage gained over manual methods by using the 
computerized system is due to the systematic pauses for 
Dranetz meter stablization prior to all measurements. 

This system has been designed to analyze any type of 
transducer. Specific adaptation to handle a particular 
transducer type with known resonance range can drastically 
reduce much of the operator interface necessary for 
evaluation of the properties of the transducer and reduce 
the time requirement accordingly. Adaptation for a 
particular type of transducer will allow variables currently 
in the computer program to be made constant. Careful 
selection cf these values will maximize accuracy and 
precision in the desired measurements and calculations. 

D. RECOMMENDATIONS FDR FURTHER IMPROVEMENTS 

In all of the preceding discussions, it must be apparent 
that the "weak link" in this system is the Dranetz impedance 
merer which is over fifteen years old. The difficulties 



39 



with this device ace primarily its lack of computer 
controllability, its plus or minus two percent accuracy, and 
its high noise at low signal levels. 

Although it was unfortunately not available for this 
study, the Hewlett-Packard i192A Low Freguency Impedance 
Analyzer seems to be the proper replacement for the Dranetz. 



It is 


capable of fall program 


control 


via 


the 


HP-IB and 


would 


also 


eliminate the neei 


for a 


sspa 


rate 


f requ ency 


synthesizer 


and data acquisiti 


on syst 


eoa. 


It 


is highly 


recommended 


that this device be 


substituted 


in 


any future 



applications of this technique. 
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APPENDIX A 



Two programs wees written to implement the measurement 
procedures and the calculations to evaluate any type of 
transducer. The first program titled "HYDRA2" automatically 
looks at a frequency range of ten times the bandwidth on 
either side of the resonance frequency. (This is modified 
for water to be five times the bandwidth on either side.) 
This program also automatically selects the frequency at 
which to take shunt measurements based on the Q factors. 
The second program, titled ’'PICXBW', allows the operator to 
specify the bandwidth for analysis and the frequency at 
which to measure shunt values. The following reference 
guide is applicable for both programs. Major differences in 
the programs occur between lines 2501 and 2650 and between 
lines 8601 and 8635. 
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REFERENCE GUIDE 



Line Nos. 


De scription 


1 - 60 


Os clarations, inputs, "bookkeeping" 


60 - 90 


Admittance Spectrum (freq. S Q) using 
subroutine 7000 


100-149 


Impedance Soectrum (freq. & Q) using 
subroutine 7000 


150 


Obtains blorked values using subroutine 
2501 


153-176 


Impedance data collection (option) 


177-186 


Admittance data collection (option) 


198-450 


Listing and plots of data using 
subroutines 


451-500 


Motional data calc, for magneto, (for 
' ? ' data) 


525-540 


Motional data calc, for elec, (for 
'!• data) 


541-642 


List of motional data and plot 


645-657 


Motional data calc, for elec, (for ' Z' 
da ta) 


658-668 


Motional data calc, for magneto, (for 
'S' data) 


670-757 


List of motional data and plot 


800-940 


Calculations for elec, coupling in air 


1000-1031 


Calculations for magn. coupling in air 


1040-1146 


Calculations for elec, coupling in water 


1 150- 1225 


Calculations for magn. coupling in water 


1350-1410 


Options for rerun or change of medium 


2501-2650 


Subroutine to determine blocked values 


7000-7545 


Subroutine for search and spectrum 
analysis 


8100-81 10 


Subroutine to print data list 


8200-8250 


Subroutine to plot real data vs. freq. 


8300-8380 


Subroutine to plot imaginary data vs. 
freq. 


8400-8450 


Subroutine to plot 'circles' 


8549-85 73 


Subroutine to obtain normalization 
fa ctors 
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8600 - 8699 
870 0-8799 
9000-9350 



Subroutine to collect data 
Subroutine to determine max/min values 
Subroutines for printouts 
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1 REM **HYCRA2** 

2 OPTION EASE 1 

3 SHORT Q0,Ol,Q8,Q9 ,Q , E ,E2 , E 3 , E5 ,R8 ,R7 ,X 2 , X 3 ,X5 , X6 , X7 , X8 
,G 7 , G8 ,£2 ,£ 3 , £5 ,£6 ,E7 ,E8 , K , Kl ,P 5 ,P6 , w 

7 SHORT A,F2 f F3,F4,F5 ,H 3 ,K 4 , J 1 , J 2 , J 3 , J4 ,R ,X,G ,E 

8 INTEGER M , Y , Y 1 , 1 , Y 2 , J , Z2 ,H 1 ,H2 ,L 2 ,S 1 

10 DIM E$ [ 2] ,A$ [2] ,TS[1] ,D (.300,3) ,A1$.[4] ,D$[30G] ,L$[20] , 
R$.[l] ,I$.[1] ,S$.[10] , CS [ 5] ,G$. [2 3] ,H$ [ 2 3] 

11 PEAL S ,N ,N8 ,N9 ,R0 ,RS ,X0 ,XS ,L ,C0 ,LG ,0 ,P1 ,H8 ,H9 ,P,W0 ,G0 
,GS ,B0 , BS ,N( 10) ,M( 10) ,F0 ,F1 ,F6 ,F8 ,F9 

14 REAL S6 ,E4,E6 ,L5 ,D ,D1 ,D2 ,D3,D4 

15 Xb ,D 1 ,D2 ,03,04=0 

16 Q,E2 ,E3 ,£ 5,R8 ,R7 ,X2,X3 ,X5=0 

17 23 ,A ,?2 ,F3 ,F4 ,F5 ,B6 ,B7 ,£6 = C 

18 S8 ,E4 ,E6 ,M,Y,Y1 ,1 ,Y2,J ,22 = 0 

19 K ,Kl,F5,P6 ,W ,00 ,Q1 ,08 ,0S ,E=0 

20 CO ,LC ,0 ,P1 ,E,KS ,K9 ,P ,W0 ,S=0 

21 N ,N8 ,N9 ,R0 ,R9 ,X0 ,X9 ,G0 ,G9 = 0 

22 80,39 ,Hl,H2 ,F0 ,F1,F6 ,F£ ,F9=0 

25 DISP "THIS PROGRAM AUTOMATICALLY SELECTS HALFWIDTH FC 
R DATA COLLECTION." 

26 DISP "THIS LOOKS AT 10 X BW IN AIR & 5 X BW IN WATER. 

II 

27 DISP "TO INPUT YOUR OWN BANDWIDTH USE "PICKEW" PRCGRA 
M INSTEAD. (HIT "CCNT')" 

29 PAUSE 

30 CLEAR § BEEP 

.35 DISP "ENTER TYPE OF TRANSDUCER TO BE MEASURED? ( 'ft ' FC 
R MAGNETIC COUPLING OR 'e" FOR ELEC.)" 

42 INPUT T§@ DISP "TO GET A COMPLETE SET OF DATA YOU NEE 
D MEASUREMENTS IN ECTH AIR AND WATER. " 

45 DISP "DC AIR FIRST. IN WHAT MEDIUM APE YOU OPERATING? ( 
ENTER ' 1 ' FOR AIR OR "2" FOR WATER)" 

46 INPUT M<s CLEAR @ BEEP 

47 IF M=2 THEN 49 

48 C $.=" AI R" ® Z2 = 10 <e CGTO 51 

49 C$=" WATER" (3 Z2 = 5 

51 PRINTER IS 2 $ PLOTTER IS 1 

55 DISP "SET THE DRANETZ FOR ADMITTANCE .SET CN LOWEST SC 
ALE THAT vvCN '.1 PEA a DURING RUN . iCGNl) " 

6 C PAUSE 

65 PRINT @ PRINT <3 PRINT 

70 PRINT "ADMITTANCE IN ";CS 

75 GGSUE 7 OCC 

80 IF M=2 THEN SO 

85 F 0= F 6 @ 00=0 <3 GCTC 100 

90 F 1= F6 <e 0 1=0 

100 GCLEAR § CLEAR @ BEEP 

105 DISP "SET DRANETZ FCR IMPEDANCE- ( SCAL E NEEDED) -HIT C 
CM " 
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lie PAUSE 

115 PRINT 0 PRINT 0 PRINT 
120 PRINT "IMPEDANCE IN ";C$ 

125 GOSUB 7000 

130 IF M=2 THEN 140 

135 F8=F6 0 Q3=C 0 GCTC 141 

140 FS=F6 0 Q9=Q 

141 CLEAR 0 EEEP 

142 PRINTER IS 701,76 0 PRINT USING 144 
14 4 IMAGE 2/ 

145 PRINT "F ( Y-AI R ) = " ; I NT (F C) , F ( Y-toAT ) = " ; INI ( F 3 ) , F 

( Z. -A I R ) = " ; IN.T ( F8) , F ( Z -NAT ) = " ; INT (F 9 ) 

146 PRINT " C ( Y-AI K) = " ; U C , w ( Y-WATEF ) = " ; Cl » * "0 ( Z-AIR) = 

";Q8 , ,"C(Z-WATER) = " ;09 

147 PRINTER IS 2 

148 CLEAR 0 DISP "THESE ARE THE INITIAL ESTIMATES. (HIT " 
CCNT' WHEN READ* TO PROCEED)" 

149 PAUSE 

150 GCSU5 2501 

152 CLEAR 0 EEEP 

153 DISP "TO COLLECT ADMITTANCE DATA (FOR ELEC. COUPLING 
) ENTER "1". ENTER "t' FOR IMPEDANCE" 

154 INPUT L20 PRINTER IS 2 

155 IF L2=2 THEN 158 

156 R $ = "G" 0 I$ = "E" 0 S$=" ADMITTANCE" <t H$= "CONDUCTANCE 

(MICRO MHOS.) " 0 G$ = " S OS CEPTANCE (MICROMHCS ) " 

157 GOTO 159 

158 RS= "R " 0 I$="X" 0 SS=" IMPEDANCE" 0 G$=" REACTANCE (OH 

MS)" 0 H $ = "RES I STANC E (OHMS)" 

159 IF L2=l THEN 176 

160 CLEAR 0 DISP "SET ON Z AND ENTER SCALE FACTCF TO COL 
LECT DATA." 

361 INPUT S60 PRINT S8;" IN ";C$ 0 PRINT "SCALE =";S8;" 
OHMS " 

162 CLEAR 0 DISP "I AM WORKING" 

164 GCSUE 6600 

165 IF Y=2 THEN 152 

166 IF M=2 THEN 169 

167 R6 = E 0 X3=E 2 0 X8=E3 0 F4=E4 

368 X6=£ 5 0 F5 = E6 0 F6 = E1 0 GCTC 198 
169 R7=E 0 X2=E2 0 X7=E3 0 F4=£4 
3 70 X 5=£ 5 0 F5=E6 0 ES=£l 0 GOTO 198 

176 CLEAR 0 BEEP 

177 DISP "SET CN Y 5. ENTER SCALE FACTOR IN MICROMHCS TO 
COLLECT DATA . " 

178 INPUT SSe PRINT S 9 ; " IN " ;C? 0 PRINT "SCALE =";S8;“ 
MICRGMhOS" 

179 GOSu'B 8600 

ISO IF T=2 THEN 152 
181 IF M=2 THEN 185 
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162 G8=E @ E3=E 2 6 E8=£3 0 F2=E4 
3 63 B6=E 5 @ F3 = E6 (§ F0=E1 @ GCTC 196 

185 G7=E @ B2=E2 & B7=E3 

186 F2=E 4 @ B 5 = E5 @ F3 = E6 § F1=E1 
156 PRINT USING 1SS 

155 IMAGE 2/ 

200 DISP "DO YOU DESIRE A LIST OF COLLECTED DATA NEAR RE 
SONANCE? (100 POINTS) ( 1 = YES ,2=NC) " 

201 INPUT Y 

202 IF Y=2 THEN 300 
204 IF L2=2 THEN 240 

210 PRINT "ADMITTANCE DATA FOR ";C$;" IN MICROMHOS" (s PR 
INT USING 215 

215 IMAGE 2/ ," FREQUENCY" , 4X , "REAL" ,6X , " IMAGINARY" 

216 GOTO 265 

240 PRINT "IMPEDANCE DATA FOR ";C$;" IN OHMS" IS PRINT US 
ING 215 

265 GCSUB 8 1U 0 

300 CLEAR Is DISP "DC YOU DESIRE A PLOT OF DATA ? ( 1= Y ES , 2= 
NO ) " 

305 INPUT Y @ GCLEAR CLEAR 
315 IF Y=2 THEN 450 

.320 DISP "ENTER: 1= G/R VS. FREQ; 2= B/X VS. FREQ; 3= B 
VS. G/X VS. R; 4= END PLOTTING LOOP." 

325 INPUT Y 
330 IF Y=2 THEN 385 
3 35 IF Y = 3 THEN 4 25 
340 IF Y=4 THEN 450 
345 GOSUE 8200 

356 LDIR 0 6 PEN 1 Is PENUP 

357 MOVE C( 1 , 1) + . 2* (D ( 30 0 , 1) -D ( 1 , 1) ) ,-(E* . 8) 

360 LABEL R$;" VS. FREQUENCY IN ";C$ 

361 GOTO 320 
365 GOSUE 8300 

357 LDIR C is PEN 1 IS PENUP 

358 MOVE D ( 1 , 1) + .2* (D(. 300,1) -D( 1,1) ) ,-L-.5*L 

400 LAEEL 15;" VS. FREQUENCY IN ";C8 

401 C-CTO 320 

425 GCLEAR Is CLEAR IS GOSUE 6400 
431 PEN 1 

442 LDIR C § PENUP 

443 MOVE - ( . S*L) , -L- . 2*L 

44 5 LABEL "INPUT ELECTRICAL. " ; S 8 ; " PLOT ( " ; 1 8 ; " vS . " ; R$; " 
) IN " ;CS 
446 GOTO 320 

450 IF L2=2 THEN 645 

451 IF T8="E" THEN 525 

h 52 GCLEAR 0 CLEAR 0 DISP "I AM CALCULATING MOTIONAL DAT 

A" 

465 C0=S* .000001/ (2*PI*£6) 
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466 LC=L/ ( (2*PI*F6) ~2*C0) 

470 FOE 1=1 TO 300 

475 D < I ,2) =D(I , 2 ).-G 

480 D (I , 3) =D (I , 2).-D (1,1) *B/F6 

485 NEXT I 

450 GOSUE S7C0 

500 GOTO 541 

525 GCLEAR @ CLEAR @ DISP "I AM CALCL'LA'I ING MOTIONAL DAT 
A" 

532 C0=B*. 00000 1/(2*FI*F) 

535 FOR 1=1 TO 300 
5 36 D (I ,2)=D(I , 2 ).-G 

537 D ( I , 3 ) =C <1,5 ).-D ( I , 1 ) *S/F 

538 NEXT I 

540 GOSUE 8700 

541 CLEAR 0 GCLEAR <s BEEP 

542 DISP "DC YOU DESIRE A LIST OF MCTICKAL DA 'I A NEAR RES 
ONANCE? ( 1=YES , 2=NC)" 

543 INPUT Y 

545 IF Y=2 THEN 604 

575 PRINT "MOTIONAL ADMITTANCE DATA IN ";C$ @ PRINT USIN 
G 580 

580 IMAGE 2/ FREQUENCY" , 4X , " REAL " , 6X ," IMAGINARY " 

595 GCLEAR £ CLEAR 
6 CO GOSUE 8100 

604 CLEAR £ SEEP @ BEEP 

605 DISP "WANT A PLOT? ( 1=YES , 2=NC ) " 

610 INPUT Y2 

615 IF Y2=2 THEN 755 
b 16 PLOTTER IS 705 0 PEN 1 
620 GOSUE 8400 

628 LDIR 0 0 PENUP 0 PEN 1 

629 MOVE - ( . 5*L) ,-L-. 2*L 

630 LAEEL "MOTIONAL " ;S$;" PLOT FOR ";C$ 

642 GOTO 755 

645 CLEAR 0 GCLEAR ® DISP "I AM CALCULATING MOIICNAL DAT 
A" 

646 IF T $ = " M " THEN 658 

647 FOR 1 = 1 TO 300 

648 D( I ,2 ) =D( I ,2) -R 

649 CG=1/(2*PI*F*X) 

650 L0=l/( (2*PI*F6)~2*C0) 

651 D (I , 3) =D (I , 3).-2*PI*D (I , 1) *L0 
b 5 2 NEXT I 

653 GOSUE 8700 

654 C0=AES (CO) 0 L0=ABS(L0) e BEEP 

655 DISP "DO YOU WANT A LIST CF THE DATA NEAR RESONANCE? 
( 1=YES ,2=NO) " 

656 INPUT Y 

657 GOTO 670 
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658 LG=X/(2*PI*F6) 

662 FOR 1 = 1 TO .200 

663 D ( I ,2 ) =D( I ,2) -R 

664 D(I ,3) =D( I , 3 ).-C ( I , 1) *2*PI*LG 

665 NEXT I 

666 GCSUB 8700 

667 BEEP @ CISP " DC YOU WANT A LIST CF DATA NEAR RESONA 
NCE? ( 1=YES ,2=NQ) " 

668 INPUT Y 

670 IF Y=2 THEN 706 

665 PRINT "MOTIONAL IMPEDANCE CATA IN ";C$ @ PRINT USING 
6 6 6 

666 IMAGE 2/ FREQUENCY" , 4X REAL ", 6X IMAGINARY " 

700 GCLEAR 5 CLEAR 

705 GCSUB 6100 

706 CLEAR @ BEEP @ BEEP 

710 DISP "WANT A PLOT? ( 1=YES ,2=NO) " 

715 INPUT Y2 

720 IF Y2=2 THEN 755 

721 PLOTTER IS 705 @ PEN 1 
725 GCSUB 8400 

737 LDIR 0 e PEN UP 

738 MOVE - (.5*L) ,-L-.2*L 

740 LABEL "MOTIONAL " ;S$?" PLOT IN ";C$ 

75 5 GCLEAR CLEAR 

756 PRINT USING 757 

757 IMAGE 3/ 

600 DISP "I AM DOING CALCULATIONS FOR TOO. PLEASE BE PAT 
I ENT . " 

603 IF M=2 THEN 1040 
oiU IF T$="M" THEN 1000 
825 IF L2=2 THEN 860 
840 w0=E 1/AES (E4-E6) 
o 45 C1=ABS (E2/ (2*PI*F0) ) 

655 GOTO 885 

860 w S=E 1/ ABS (E4-E6) 

db 5 C 1=A£S ( 1/ ( 2* P I * F8 *E 2 ) ) 

685 K= 1- (F0/F6 ) "2 
695 Kl=Cl/(Afc£ (CO) +C1) 

33G GCSUB 9000 
921 Pf INTER IS 2 
940 GOTO 1350 
1000 K=l- (F6/F0) "2 
1005 IF L2=2 THEN lo25 
1010 U 0=E 1/AES (E4-E6 ) 

1015 GOTO 1 U 30 

1025 <*.£ = £ 1/AES (E 4-E6 ) 

lO 30 GCbUE 92 Go 

1 U 3 1 PRINTER IS 2 e GOTO 1350 

Id 40 IF T$> = "M" TEEN 1150 
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1041 DISP "If YOU HAVE DATA FOR AIR, ENTER 1; ELSE 2" 

1042 INPUT Y 

1043 RAD 

1044 IF L2= 1 THEN 1110 

1045 IF Y=2 THEN 1062 

1060 P6=D4* (D3-D4)/(D3*R7) 

1061 vv0=Fl 

1062 Q1=E 1/ABS (E 4-E 6) @ K = 1- ( F 1/F9) "2 
1070 GOSUB S100 

1076 PRINTER IS 2 <s GOTO 1350 
1110 IF Y=2 THEN 1135 
1126 WC=F 1 

1130 P6=D2’ f (D1-D2)/ (Dl*G7) 

1135 C1=E1/AES (E4-E6) 

1140 K=l- (Fl/FS) *2 

1145 GCSUB 5100 

1146 PRINTER IS 2 @ GGTC 1350 

1150 RAD 

1151 DISP "IF YOU HAVE DATA FOR AIR, ENTER 1; ELSE 2" 

1152 INPUT Y 

1153 IF L2= 1 THEN 1160 

1154 IF Y=2 THEN 1162 

1155 B=.5*ACS( (R7-R9.)/D4) 

1156 H b=S IN (E) "2 g H 9=COS (£) ~2 

115 7 P 5= ( ( 1+D 3/R0*H 5) 5- ( L-D 3/R0*H8 ) ~.5)/( (1+D3/RC*H9) 

. 5+ ( 1-D 3/R0 *H8 ) * .5) 

1155 F=D4*SIN (2*d)/(4*R9*CS) @ wG= F 9* ( 2*P + ( P "2+ 1) 3 . 5) 
1160 P6=L4 * (D3-D4 ) / (D 3*R7 ) 

1162 C9=E 1/ABS (E 4-E 6) Q K= l-( F 5/F 1 ) * 2 

116 3 PRINT "£= " ; E , " K= " ; K , " Kl= " ; Kl , " K£ = " ; H8 ," H 9= " ; H5 
1164 GCSUE 9300 

1168 PRINTER IS 2 6 GOTO 1350 
1160 RAD 

1162 IF Y=2 THEN 1210 
1184 £=. 5*ACS ( (G7-G9 )/D2) 

1186 H8 = SIN (B) A 2 @ H9=CCS (B) "2 

1186 P 5= ( ( 1+D 1/G0*H 9 ) 5-( 1-D1/G0*H8) 5)/ ( ( l+Dl/Gu*H 5) 

. 5+ ( 1-D 1/G0 *H&) * .5} 

1190 P=D2*SIN ( 2* B) / ( 4 *G9*v 1) Is *0=F1* (2*P + (F "2+1) .5) 

1 lb? 5 Pb=D2* (D1-D2)/ (D1*G7) 

1210 C5=E 1/AES (E4-E6) 0 K= 1-(F 9/F 1) ~2 

1211 PRINT " E = " ; E ,'"K= " ; K ," K 1=" ; K 1 , " H8 = " ;Hb ,"H5=";H9 
1215 GOSUB 5300 

1225 PRINTER IS 2 @ CLEAR is ELEP 

1350 DISP "DC YOU DESIRE ANOTHER IN THIS NED IUM? (TO GET 
THE OTHER TYPE DATA? ) ( 1=YES , 2=NO ) " 

1355 INPUT Z1 
1357 PLOTTER IS 1 

1360 IF Zl=l THEN 152 

1361 PRINT USING 1362 
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1362 IMAGE 3/ 

1365 DISP "DC YCU DESIRE A RUN IN ANOTHER MEDIUM? ( 1=YES , 
2=NG ) " 

1370 INPUT Z1 

1375 IF Zl=l THEN 45 

1380 DISP "JUST FOR THE RECORD' , INPUT WATER TEMP, AIR TE 
MP, TRANSDUCER SER NO, MODEL NO.," 

1381 DISP "AND TYPE (E OR M ) " 

L 28 5 DISP "INPUT *'C" IF INFO IS UNKNOWN" 

13 50 INPUT W ,A ,S ,N ,L6 

1 3 S 1 PRINTER IS 701 ,76 

1400 PRINT "WATER TEMP=" ; w , ,"AIR TEMP= " ; A , SER . NG. = ";S 
, ," MODEL NO . = " ;N , ," TYPE " ; L$ 

14 10 GOTO 5958 
2501 REM *ShUNT* 

250 2 CLEAR § BEEP i§ DISP "WE ARE FINDING SHUNT VALUES FG 
R G/E OR R/X." 

2510 A $ = " F R " <e ES.="HZ" @ Sl=l 

2511 IF T$="M" THEN 252G 

2512 IF M=2 THEN 2516 

2512 Q=MAX(Q0,Q5) @ F=FG- (Z 2-2) *F0/Q 

2514 IF F>2C0 THEN 2535 

2515 F-20G+ (FG-2*FG/Q-2CG)/5G 0 GOTO 2535 

2516 Q=MAX (Q1,Q9.) (° F=F ]- (Z 2-2 ) *F 1/y 

2517 IF F>200 THEN 2535 

251c F=20U+(F1-2*F1/Q-2U0)/5G § GOTO 2535 

2520 IF M=2 THEN 2524 

2521 Q=MAX(Q0,Q8) 0 F=Fd- (Z 2-2) *F8/Q 

2522 IF F>200 THEN 2535 

2523 F = 20C+ ( F8.-2* Fb/Q-2 00) / 50 0 GOTO 2535 

2524 0=MAX (v 1 , t 9.) 0 f =F S- (Z 2-2) *FS/Q 

2525 IF F>200 THEN 2535 

2527 F=20G+ (F9-2*F9/Q-2CG)/ 50 
2525 D$=VAL$ (F ) 

2527 OUTPUT 717 ;A S,DS,£$ 

2538 WAIT 1000 

2540 DISP "SET DRANETZ FREQ. SCALE TO COVER FP.Ev . ON SYN 
THESIZER .ZERO METERS .SET ON FS.(CCNT)" 

2541 PAuSE 

25«2 DISP "HEPE aE GET THE NORMALIZATION FACTORS." 

2543 OUTPUT 709 ;"vTl" 

2546 GGSUB b54S 

2550 G,B,R,X=0 

2551 A 8 = "FR" 0 E$="HZ" 0 D$=vALS(F) 

2553 OUTPUT 717 ;A$,D$,B$ 

2561 CLEAR <5 DISP "SET DRANETZ ON Y ,SET SCALE , SET ON NOR 
>1 FGF FILTER + PHASE, PLUG IN TRANSDUCER." 

2562 BEEF g DISP "ENTER SCALE FACTOR IN M I CR CM EG S . S Am E S 
CALE AS FOR DATA . (v«E GET SHUNT G/B)" 

2564 INPUT Set? CLEAR 0 DISP "I AM aGRKING TO GET S HUNT v 
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ALoES" £ beep 

2565 FOR 1=1 TO 10 

2566 OUTPUT 717 ?A§ ,D$,E$ 

2567 WAIT 1000 

2 566 OUTPUT 705 ; "AC 3 " 

2569 ENTER 705 ; M ( X ) 

2570 G=G+K(L) 

2571 OUTPUT 709 ;"AC4" 

2572 ENTER 709 ; N(I) 

2573 B=B+N ( I) 

2574 NEXT I 

2575 G=G/ ( I — 1 ) £ E=B/(I-j) 

2576 G=G*S3/K6 £ B=B*S6/NS £ BEEP 

2576 CISP "SET CRANETZ ON Z ; SET SCALE FOR LAX RESPONSE. fc 
CVE TRANS. INPUT. ENTER SCALE FACTOR " 

2575 INPUT ScC« CLEAR £ DISP "I AN NO F KING TO GET SFUNT 0 
ALOES" 

2 560 FOR 1=1 TO 10 

2581 OUTPUT 717 ;A$,D$,E$ 

2582 aA IT 1000 

2563 OUTPUT 709 ;"AC3" 

2564 ENTEF 70S ; N(I) 

256 5 R=R+i*i(L) 

25c 6 CUTFL.T 709 ;"AC4" 

2567 ENTER 705 ; N(I) 

2566 X=N ( I ) +X 
2585 NEXT I 

2550 R=R/(I-1) £ X=X/(I-j) 

2591 R=R*S8/N8 £ X=X*S6/N5 

2592 IF T $= "M" THEN 2610 
25 9 5 IF N= 1 THEN 26 Cl 

2556 R5=R £ XS = X £ G9 = C £ BS=B £ GOTO 2601 
2 6ul G0=C 6 E0 = E & R0=R £ Xu = X 

2601 PRINT "SHUNT VALUES " 

2602 PRINT USING 2603 



2 60 ; INAGE 1/ 

2605 PRINT "G0 = " ; GO , , " £0 = " ; £ 0 , ," P0= " ; RG 
;GS, L 5= " ; B9 , ,"R9=";R9 , ,"X9=";X9 



," X0 = " ; XC , 



u 



2 6 06 PRINT is GCTC 26 30 
2o 10 IF S 1=2 THEN 2620 

2611 J 1=G £ 02= £ (s 0 3=R y J4 = X 

2612 IF N = 2 THEN 2ol5 

2613 F =F6+ ( Z2-2 ) *F 6/v § Sl=2 £ GCTC 2550 
2615 F=F9+ (Z2-2) *F9/U £ Sl=2 £ GCIC 2550 

2620 G= . 5* (G+0 ] ) is B=.5*(S+J2) £ R=.5*(F+U3) y X=.5*(U4+ 

X) £ GCTC 2595 

2630 GCLEAR £ CLEAR y E EE P 

2 c 3 5 CISP "ENTER ' V IF ALL IS WELL; "2' IF ZCl NEED A R 
EPEAT 

2640 INPUT \ 
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2645 IF Y=2 THEN 2501 
2650 RETURN 

7000 REM **F-Q** 

7001 GCLEAR @ CLEAR 

7005 D ISP 11 ENTER LOWER, UPPER FREO . FOR SWEEP. IF PESCNAK 
CE UNKNOWN USE 10000,200000. (Fl,F2)" 

7010 INPUT F 4 , F 2 

7011 DISP "ENTER RMS VOLTAGE NEEDED. ( 3000 MV FOR LOW FRE 
Q . AS LOW AS 2000 MV FOR HIGH)" 

7012 INPUT A§ A 1$ =VAL$ (A ) 

7020 CLEAR <f BEEP 

7025 A $ = "AM" 3 3$ = " MR" 

7035 OUTPUT 717 ;A?,A1$,B$ 

7040 OUTPUT 709 ;"AC3VT3" 

7045 F3=CEIL( (F2-F4) /3CC) 

7060 CLEAR <g DISP "I AM WORKING TO GET THE SPECTRUM DATA 
FOR YOU , 1HEN WILL SHOW & MAKE A PLOT" 

7065 FOR 1=1 TO 300 
7070 D ( I , 1 ) =F4 + I *F 3 
7075 A $= "HZ" £ B$. = "FR" 

7063 DS=VAL$ ( D ( I ,1) ) 

7085 OUTPUT 717 ;B$,D$,A$ 

7090 OUTPUT 709 ;"AC3VT3" 

7095 WAIT 100 

7100 ENTER 705 ; D(I,2) 

7105 NEXT I 

7110 F5=F4-. 1* (F2-F4) <s F 6=F 2+ . 1* (F 2-F 4 ) 

7115 F6=F2+» .1* (F2-F4 ) 

7119 PLOTTER IS 1 § G CL EAR <s CLEAR 

7125 SCALE F5,F6 ,-.l ,1.2 

7130 XAXIS 0,2000,F4,F2 

7135 YAXIS F 4 , .1 ,-.l ,1.2 

7140 FOR 1=1 TO 300 

7145 PENUP 

7150 PLOl D (I , 1) ,D( I ,2) 

7155 NEXT I 
7158 GRAPH if COPY 

7160 DISP "NEED ANOTHER FUN? ( 1= YE S , 2=NC ) " 

7161 INPUT * 

7162 IF 7=1 THEN 7005 

7163 CLEAR fe GCLEAR vf BEEP 

7165 DISP "ENTER DECISION POINT FOR AMPLITUDE (0 TO 1.2) 
. (LESS THAN THE MAX DISPLAYED) ." 

7170 INPUT Sts CLEAR $ BEEF 

7180 PFINT “AMPLITUDE IN VOLTS";" FREQUENCY" 

7161 PRINT 

719 0 FOR 1 = 1 TO 20 0 

715 5 IF C ( I , 2 ) <S THEN 7210 

72 UU PRINT USING 7 20 5 ; D ( I , 2 ) ,D ( I , 1) 

7205 IMAGE IX , D . DDDDDDC D , 1 OX .CDDDDD.LD 
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7210 NEXT I 

7215 PRINT USING 7220 

7220 IMAGE 3/ 

7225 CLEAR (e BEEP 

7226 DISP "IF YOU ARE READY TC CONTINUE ENTER 'l'; OTHER 
WISE ENTER ' ' t' 

7227 INPUT Y 

7228 IF Y=2 THEN 7160 

7235 DISP "WHAT IS CENTER FREQ, HALFWIDTH TC CONSIDER?" 

7236 INPUT F,N@ CLEAR 0 BEEP 

7240 DISP "I AM WORKING TO GET GOOD FREQ AND Q DATA" 

7241 FOR 1=1 TO 300 

7242 D ( I , 1 ) = F. -K+ 1 * N/ 1 5 0 

7243 A $ = " H Z " 0 B $ = ” F R ” 

7245 D $=VAL$ ( D 11,1)) 

7255 OUTPUT 717 ;B$,D$,A$ 

726 C OUTPUT 709 ;"AC3VT3" 

7265 WAIT 100 

7270 ENTER 70S ; D(I,2) 

7275 NEXT I 

7285 CLEAR 0 GCLEAR 0 DISP "I AM FINDING THE ABSOLUTE MA 
X AND FREQ UPPER AND LOWER" 

7 2S0 £1,B4 = 0 0 H ,H 1 ,H 2 ,H 3 , H 4= 1 

72S5 E2,£3=5U 

7300 FOR 1=2 TC 300 

7 310 IF C (I ,2) <D (H ,2) THEN 7 340 

7315 IF D (I ,2) =D (H ,2) THEN 7 325 

7325 A6 = D (1,2) 0 F6 = D(I,1) 0 H=I 

7340 NEXT I 

7345 A7=A6/2 

7350 FOR 1=1 TC H 

7355 IF A7=D( I ,2) THEN 7405 

7360 IF A7 <D ( I ,2) THEN 7 365 

7 36 5 IF E ( I , 2 ) <E 1 THEN 7 42 5 

7 370 B 1=D (1,2) 0 H 1 = 1 0 GOTO 7425 

7365 IF D ( I , 2) >B2 THEN 7400 

7 3 SO B2 = D ( I ,2) 0 H2 = I 

7400 GOTO 7425 

7405 E1,B2=D<I,2) 0 Kl,H2=I 0 F7=D(H2,3) 0 GOTO 7440 
7425 NEXT I 

7430 X= (A7-B])/(B2-B1) 

74 3 5 F7=X* ( D ( H 2 , 1 ) — C (Hi ,1))+C ( ti 1 , 1 ) 

7440 FCR I =H TC 300 
7445 IF A7=C (1,2) THEN 74S5 
7450 IF A7>D( I ,2) THEN 747 5 
7 455 IF D ( I ,2) >B3 THEN 7 47 0 
7460 B3=D (I ,2) 0 H 3= I 

7470 GOTO 7515 

7475 IF D ( I ,2) <=£4 THEN 7490 
7480 B4=C ( I ,2) 0 H4=I 
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7490 GOTO 7515 

7495 B3,34 = D(I,2) @ H2,H4=I 

7505 L5=D(H2 , 1) 

7510 GOTO 7530 
7515 NEXT I 

7520 X=(A7-B4)/(E2-B4) 

7 525 L5=- (X* (D(H4 ,1).-D (H 3 , 1 ).) ) +D ( H4 , 1 ) 

7530 Q=F 6/ ( L5-F7) 

7535 PRINT "CENTER FREQ IS " ;F6 

7536 PRINT "C IS ";G 

7538 BISP "IF YOU ARE READY TO PROCEED, ENTER ' 1 ' . TO RE 
RUN FOR EET'TER VALUES, ENTER "2'." 

7539 INPUT Y 

7540 IF Y=2 THEN 7235 
7545 RETURN 

8100 REM **DATA LIST** 

8101 PRINTER IS 2 

8102 FOR 1=99 TO 199 

8103 PRINT USING 8105 ; D (I ,1) ,D(I,2) ,D(I,3) 

8105 IMAGE DDDDDD ,D,2X ,D.DCDE ,2X,D.DDDE 

8106 NEXT I 

8107 PRINT USING 8106 
810b IMAGE 3/ 

6110 RETURN 

8200 REM * *PLCT RE** 

8201 GC LEAR <s CLEAR 

8202 H8*D( 1,1)-. 15* (0(300 ,1)-D( 1,1) ) 

620 3 H9=D (300 , 1) + . 1* (D ( 300 , 1) -C ( 1 , 1 ) ) 

8205 PLOTTER IS 705 <s PEN 1 

8209 SCALE H8 ,H 9 . 8*E ) , E + . 2 5*E 

6210 XAXIS 0 , (D( 300 ,1 )-D ( 1, 1) )/10 ,D ( 1 , 1) ,D (.200 ,1) 

8211 YAXIS D( 1,1) ,E/5 ,-(. 8*E ) ,E + .2*E 

8212 PENUP 

8215 FOR 1=1 TO 300 

8216 PLOT D(I , 1) ,C( I ,2) 

8217 NEXT I 

8216 PENUP @ DEG @ LDIR C,SIN(90) 

8221 FOR LL=D ( 1 , 1) TO 0(200,1) STEP ( D ( 200 , 1 ) -D ( 1 , 1 ) ) / 10 
!§ PENUP 

8222 MOVE LI,- (. 18*£) 

8 22 2 LABEL IN.T(Ll) 

8224 NEXT LI 

8226 LDIR 0 is PENUP 

8227 MOVE D ( 1, 1) + .2* (D ( 200 ,1) -O ( 1 ,])),-(. 7*E) 

6228 LABEL "FREQUENCY (HE) " 

6 2 20 PENUP is PEN 1 
6222 LDIR 0 

8224 FOR LL=- ( . 6 *E ) TO E+.2 *£ STEP E/5 £ PENUP 
82 2 5 MOVE D ( 1 , 1) -.09* (D ( 2U0 ,1) -D ( 1 ,1) ) ,L 1 
6 2 26 LABEL IN.T(L1) 
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6 237 NEXT LI 
6240 LDIR 0 ,SIN ( 90) 

3 24 1 MOVE D ( 1 , 1 ) - .1* (D (300 , 1 ) -D ( 1,1) ) 1*E 

8243 LABEL H$ 

6250 RETURN 
8 300 REN * * P LOT IN** 

6301 GCLEAR 0 CLEAR 
8 30 3 H6=D( 1 15* (0 ( 30 0 ,1) -D ( 1 ,1) ) 

6 304 H5=D ( 30 0,1) + .1*(D (300,1) -D(l,l)) 

8 305 L=MAX (AES ( E5 ) ,AtS(E3)) 

8306 PLOTTER IS 705 0 FEN 1 

6312 SCALE H8 ,K9 ,-L- . 5*L ,L+ . 25*L 

6313 XAXIS 0 ,(D( 300 ,1 )-D ( 1 , 1) )/lU ,D( ] , 3) ,D ( 300 ,1) 

8314 YAXIS C(1,L) ,L/5 ,-L - . 2*L , L+ . 2*L 

8315 PENUP 

8318 FOR 1=1 TO 300 
8 325 PLOT C (I ,1) ,D( I , 3) 

8330 NEXT I 

8 3 31 PENUP 0 DEG 0 LDIR C,SIN(90) 

83 33 FOR L1=E (1,1) TO D( 300 , 3) STEP (D ( 30 0 , 1 ) -D ( 1 , 1 ) ) / 10 
0 FEN UP 

8 335 MOVE Ll ,-L-.2*L 
6 3 36 LABEL IN.T(L1) 

6337 NEXT Ll 

8340 LDIR 0 0 PENUP 

8 34 1 MOVE D ( 1 , 1) + .2* (D(30G,1)-D (1,1) ) ,-L-. 4*L 
8342 LABEL "FREQUENCY (HE)" 

8345 PENUP 0 PEN 1 0 LDIR 0 

6352 FOR LL=-L-.2*L TO L+ ,2*L STEP L/ £ 0 PENUP 

8355 MOVE D ( 1 , 1 ) - . 09 * (D ( 300 , 1 ) -D ( 1 , 1 ) ) ,L 1 

8356 LABEL INT(Ll) 

8357 NEXT Ll 

8360 PENUP 0 PEN 1 0 LDIR 0,SIN(90) 

8 361 MOVE D(1,1)-.1*(D (300,1) -D (1,1) ) ,- (L* . 5 ) 

6368 LABEL G$ 

6380 RETURN 
6 400 REM * *C IRCLES * * 

6401 GCLEAR 0 CLEAR 

8402 L 1=MA X ( ABS ( E 3 ) ,A ES ( E 5 ) ) 

6403 L =MA X (AES (E) ,L1) 

6404 FLOTTER IS 705 0 PEN 1 

8406 SCALE 1 . 25* (-L- . 1*L ) , 1 . 25* (L+ . ] *L) ,-L- . 2*L ,L+ . 15*L 
6407 XAXIS 0 ,L/5 ,-L-.l*L ,L+. 1*L 
8408 YAXIS 0 ,L/5 ,-L-. 1*L ,L+ . 1*L 
6409 PENUP 0 PEN 1 

8411 FOR 1=1 TO 300 

8412 PLOT D (I ,2) ,D( I , 3) 

8413 NEXT I 

8414 PENUP 

6415 LDIR 0 0 PEN 1 0 PENUP 
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8 416 MOVE - (L*. 5.) ,-L-. 1*L 
8416 LAE EL H$ 

8420 PEN UP @ PEN 1 @ DEG 

8425 LDIR C,SIN(9CL) 

8426 FOR L1=-L-.1*L TO L+.1*L STEP L/ 5 (9 PENUP 
6427 MOVE LI ,-L 

8428 LAEEL IN.T(L1) 

8429 NEXT Ll 

8435 LDIR 0 @ PEN 1 @ PENUP 

8436 FOR L1=-L-.1*L TO L+.1*L STEP L/S <s PENUP 

6437 MOVE - ( .6S6*L) , Ll 

6438 LAEEL IN.T(L1) 

6439 NEXT Ll 

8445 LDIR C,SIN(SG) 

64 + 6 MOVE - ( . S5*L) ( . 3*L) 

6448 LAEEL G$ 

8450 RETURN 
8549 REM **NCRM** 

6550 N8 ,N9=0 

6551 FOR 1=1 TC 10 



6557 WAIT 5U 0 

8558 OUTPUT 709 ;"AC3" 

8559 ENTER 70S ; N(I) 

8560 N8=N8+N(L) 

8565 OUTPUT 709 ; ,, AC4'’ 
6566 ENTER 709 ; M(I) 
8567 N S = N 9 +M ( L ) 

6568 NEXT I 
8 5bS N 8 = N 6 / ( I - 1 ) 

8570 NS=-(KS/(I-1) ) 
8572 PRINT USING 6573 
8 57 3 IMAGE 3/ 

8575 RETURN 



6600 REM ** VALUES * * 



6601 CLEAR @ 

6602 DISP "I 

6605 IF L2= 1 

6606 IF L2= 1 
6507 IF L2=2 
86 Oo IF L2=2 
6620 Fo=F 1 $ 
8b 2 1 F o=F0 Q 
86 22 F 6=F6 0 
86 2 3 Fb=F 9 (s 
8625 FCR 1=1 



EEEP 

AM CCLLECTING THE DA I A NGu" 

AMD M=2 THEN 6620 

AND M=1 THEN 8621 

AND M=1 THEN 6622 

AND M=2 THEN o622 

y=vl (9 GCTC bo2 5 

O=w0 & GOTO 6625 

W=V6 (9 GOTO ti625 

C=u9 

TC 50 



6626 IF F 6 — Z 2* b 6/0 + ( Z 2-2 ) * F6 / ( 5 0 *^ ) > 2 0 0 THEN 
do 29 D (I ,1 ) =200+ (Fb-2*F6/w -200) *1/50 
6 6 30 GOTO do 32 

86 3 1 C ( I , 1 ) =F b-22*F 6/0.+ I * ( 2 2—2 ) *F6/ (50*0) 
6b32 D8=vALd (D ( I , 1 ) ) 



o b 3 1 
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8633 OUTPUT 717 ;A$,D$,E$ 

8634 WAIT 150 

86 3 5 OUTPUT 709 ; ,, AC3VT3" 

8636 ENTER 709 ; H8 

8637 C(I ,2) =H8*S8/N8 

86 38 OUTPUT 709 ;"AC4VT3” 
8639 ENTER 709 ; H5 
06 40 D ( I , 3 ) =H 9*S8/NS 
8641 NEXT I 



6650 FOR 1=51 TO 250 

86 51 C(I ,1) =F6-2*F6/Q+ ( I - 50 ) * 2* F6/ (10 0*0) 

8652 D$=VAL$ ( D ( I ,1) ) 

8653 OUTPUT 717 ?AS.D$,B$ 

8654 WAIT 250 

8655 OUTPUT 709 ;"AC3VT3" 

6656 ENTER 709 ; H8 

8657 D ( I ,2) =HS*S8/N8 

8658 OUTPUT 709 ;"AC4VT3" 

86 59 ENTER 70S ; K9 

8660 D (I ,3) =H9*S8/N9 

6661 NEXT I 

8670 FOR 1=251 TO 300 

86 71 D(I,l)=F6+2*F6/'Q+(I-2 50) * (Z 2-2) *F6/ ( 50*0) 
8672 D $=v AL$ ( D ( I ,1) ) 

6673 OUTPUT' 717 ;A 

8674 WAIT 150 

8675 OUTPUT 709 ;"AC2VT3" 

8676 ENTER 7US ; H8 
86 77 D (I ,2) =H8*SS/N8 

66 78 OUTPUT 709 ;”AC4VT3" 

8679 ENTER 70S ; H9 
o68C D (I , 3) =H9*S8/N9 

8681 NEaT I 

8682 CLEAR g DISP "I AN FINDING PiIN/NAX VALUES 
FREQS" 

8663 C-CSUB 6700 

6684 IF L2=2 THEN 8687 



6685 PRINT " GN AX = " ; E , , " 
f " FEN;X=" ;£« , uf.I=‘ ;E5 , 
8686 GOTO 8b 97 



FGMX= " ; £ i , , 
f "FBM = ";E6 



£GKX=" 






8687 PRINT "RnAX=";E, FRNX= " ; E 1 , XFF:X=" ; E2 , , 
FXMX=" ;E 4 , XMI = " ; £5 , ," FXMI. = " ; E6 
8697 DISP "IF ALL IS wELL , ENTER "l". ENTER " '2 
E DATA." 



869o INPUT Y 
6699 RETURN 



8 700 REN * *FAXNIN * * 

8 70 3 h ,H 1 , H2= 100 

8 705 E=D(lUO,2) £ El ,E4 ,£6 = D ( 100 ,1) E2,E3,£5 

870o FOR 1=101 TO 20 u 



& ASSOC. 

1 ENA = " ; E 3 , 

’ XNX = " ; E 3 , 
' TO RETAK 

=d (lot, :-) 
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8710 IF C ( I , 2 ) <0 (H 1 2 ) THEN 6725 

6 711 IF D (I ,2) =E(H ,2) THEN 87 2 0 

8720 E=D (1,2) @ El=D(I,l) @ E2=D(I,3) § H = I 

8 7 25 IF D (I ,3) <D(H1 ,3) THEN 67 40 

8726 IF D (I ,3) =D (Hi ,2) THEN 6735 

6735 E3=D ( I , 3 ) @ E4=D(I,1) g H 1=1 

6740 IF D (I ,3) >=D(H2 r 3) THEN 6751 

8745 H 2=1 @ E5=D(I,3) § E6=D(I,1) 

6751 NEXT I 

8760 J 1=D (H , 2 ) <s J2=(D (H+1,2) -0 (H-l ,2) )/2 0 J3=(D(H + 1,2) 

+D(H-l,2).-2*ul)/2 

8762 F6=-(J2/(2*J3) ) 

6 7o 4 A=J 1+02+F6+J 3*F 6 2 



o 770 F 6=D (H ,1)+F6* (D (H ,1) -C (K-l, 1) ) 

877 5 £ 1=F 6 fs E=A 
8788 D=ABS ( E 3-E 5) 

8 76 S IF k=2 THEN 87 5 3 
6750 IF L2= 1 THEN 8752 
8791 C3=D § GOTO 6756 
8 792 D 1=C e C-OTO 8796 
6793 IF L2= 1 THEN 8795 
6754 D 4= C @ GOTO 6756 
8795 D2=D 

6756 PRINT "Dl=" ;D1, ,"D2=";D2 , ,"B3 = ";D3 , f "D4 = " ;D 4 

8759 RETURN 

3000 REM **CALC-A** 



5001 PRINTER IS 701,76 

5002 PRINT g PRINT 

9003 PRINT "VALUES MEASURED IN AIR FROM " ;S5;" DATA" 
SGu4 PRINT 

SOU PRINT "BLOCKED CAPACITANCE- " ; A£S (Co ) ; " FARADS AT 
; F ; "HZ" 



9012 PRINT 

9013 PRINT 
5014 PRINT 
9015 PRINT 
SOlo PRINT 



"ELECTRICAL vUALITY FACTOR = ";G0 
"MECHANICAL yUALITY FACTOR = ";w6 
"ELECTRICAL RESONANCE =";F0;" HZ" 

"MECHANICAL RESONANCE = ";F8;" HZ" 

"DYNAMIC ELECTROMECHANICAL COUPLING COEFFICIE 



NT = 
5017 
9021 
50 31 
50 3 ^ 
9 0 £u 



" ; K 

PRINT "STATIC COUPLING COEFFICIENT - " ;Rl 
PRINT "BLOCKED RE'S I STANCE = " ; PC; "CHI- S" 
PRINT USING 9032 
IMAGE 2/ 

RETURN 



910 0 REM **CALC-vv** 

5101 PRINTER IS 701,76 
9102 PRINT g PRINT 

9 10: PRINT "VALUES FOR i- ZASU RELENTS IN 4 ATE R FROM ";S5; 
" DATA" 

9 lu 4 PRINT 

5111 PRINT "ELOCKEL CAPACITANCE = ";A£S(C0);" FARADS AT 
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9112 PRINT 
Sill PRINT 

5114 PRINT 

5115 PRINT 

5116 PRINT 
NT = " ; K 
5120 PRINT 



"ELECTRICAL QUALITY FACTCR = " ; Q ] 

"MECHANICAL QUALITY FACTOR = ";QS 
"ELECTRICAL RESONANCE = ";F1;" HZ" 

"MECHANICAL RESONANCE = ";F9;" KZ" 

"DYNAMIC ELECTRCMECHANICAL COUPLING COEFFICIE 

"MECHANICAL POKER UTILIZATION FACTOR" ; AES ( P6 ) 



5121 PRINT " STATIC COUPLING COEFFICIENT*" ;K 1 
9122 PRINT "ELCCRED RESISTANCE = " ; R5; "OHMS " 
5132 PRINT USING 9133 



9133 IMAGE 3/ 

9150 RETURN 

S2C0 REM **CALC-A/M** 

9201 PRINTER IS 701,76 

9202 PRINT (3 PRINT 

9203 PRINT "VALUES MEASURED IN AIR FROM " ; S 5 ; " LATA 



9204 PRINT 

9211 PRINT "BLOCKED INDUCTANCE- " ;AES(LC);" HENRIES AT 
;F6; "HZ" 

9212 PRINT "ELECTRICAL QUALITY FACTCR = ";o0 

9213 PRINT "MECHANICAL QUALITY FACTCR = " ; *S 



II 



9214 PRINT 
5215 PRINT 
9216 FRINY 
NT = " ; K 
5217 PRINT 



"ELECTRICAL RESONANCE =";F0;" HZ" 

"MECHANICAL RESONANCE = ";F8;" HZ" 

"DYNAMIC ELECTRCMECHANICAL COUPLING COEFFICIE 



" E LOCKED RESISTANCE = 



"CHMS 



9231 PRINT USING 9232 



5232 IMAGE 2/ 

9250 RETURN 

5300 REM **CALC-K/M** 

9301 PRINTER IS 701,76 

9302 PRINT £ PRINT 

5 303 PRINT "VALUES FCR MEASUREMENTS IN vmATER FRCM ";S$; 
" DATA" 



9304 PRINT 
£311 PRINT 
";Fo;"hZ" 

? 312 PRINT 

9313 PRINT 

9314 PRINT 
5315 PRINT 
9316 PRINT 
NT = " ; K 

5 z 1 o P R I NT 
9319 PRINT 



"£ LOCKED INDUCTANCE = ";AES(L0);" HENRIES AT 

"ELECTRICAL QUALITY FACTCR = " ; Q 1 
"MECHANICAL vUALITY FACTOR = ";Q£ 

"ELECTRICAL RESONANCE = ";F1;" HZ" 

"MECHANICAL RESONANCE = ";F9;" HZ" 

"DYNAMIC ELECTRCMECHANICAL COUPLING COEFFICIE 

"POTENTIAL EFFICIENCY = " ; t 5 

"FREQUENCY CF OPTIMUM EFFICIENCY = ";.vo;" HZ" 



5320 PRINT 



"MECHANICAL PC.vER UTILISATION 



FACTCR" ; AcS (Pt>) 
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(N| CO 



PFIN'I USING 9 32 3 
IbAGE 2/ 

9350 PET URN 

9996 DISP " THE END" 

S99S END 



110 



1 REM **PICK5W** 

2 OPTION BASE 1 

3 SHORT QG ,01 ,08 ,09 ,Q ,E ,E2 ,£ 3 ,E5 ,RS ,R7 ,X2,X3 ,X5 ,X6 ,X7 ,X8 
,G7,G8 ,E2 ,B3,B5,B6 ,B7 ,B8 , R , K 1 ,P 5 ,P6 ,v\ 

7 SHOFT A ,F2,F3 ,F4 ,F5 ,H3 ,H4,J1 ,J2 ,J 3, J4 ,G,E ,R ,X 

8 INTEGER M , Y , Y 1 , 1 , Y2 , J , 2,2 , 2 3 , Hi ,H2 ,L 2 ,S 1 

10 DIM B$ [ 2] ,A$[2] r T$[ L] ,D (300, 3) ,A1$.[4] ,D$[30C] ,L$.(20] , 
R5.[l] ,I$[1] , S $ [ 10 ] ,C$[5] ,G$ [2 3] , H $. [ 2 3 ] 

11 REAL S ,N ,N8 ,N9 ,RG ,R9 ,X0 ,X9 ,L ,C0 ,L0 ,0 ,P1 ,H8 ,H9 ,P ,WG ,G0 
,G9 ,E0 ,B9 ,N ( 10) ,M(1Q) ,F0,F1,F6,F8,FS 

14 REAL S3 , E4 , E6 , L 5 ,D ,D1,C2 ,D 3 , D4 

15 X6,Dl r D2,D3,C4 = 0 

16 C,E2 ,E3 ,E5 ,R8 ,R7 ,X2,X3 ,X5=0 

17 Z 3 ,A ,F2 ,F3 ,F4 ,F5 ,B6,E7,B8=C 
16 S8,E4,E6 ,M , Y , Y 1 , 1 ,Y 2 , J , Z 2=0 

19 K ,K1 ,P5 ,P6 ,W ,00 ,Q1 ,08 ,09 ,E = 0 

20 CO ,LG ,C ,P1 ,B,H6 ,H9 ,P ,W0 ,S=0 

21 N ,N8,N9 ,R0 ,R9 ,X0 ,X9 ,G0 ,G9=0 

22 BO , ES ,H 1 ,H2 ,F0 ,F1 ,F6 ,F8 ,FS=0 

24 DISP ".THIS PROGRAM IS DESIGNED FOR AN EXPERIENCED CPE 
RATCR. THE BANDWIDTH FOR DATA COLLECT." 

25 DISP "AND FREC. TO GET E LOCKED DATA WILL EE THE CHCIC 
E OF THE OPERATOR. "HYDRA 2' PRCGRAM IS" 

26 DISP "DESIGNED TC CCVER 10 X BW AUTOMATICALLY. HIT ' 'C 
OUT" IF READY TO FRCCEED . " 

27 PAUSE 

28 CLEAR @ BEEP 

.35 DISP "ENTER TYPE OF TRANSDUCER TC BE MEASURED? ( 'M ' FC 
R MAGNETIC COUPLING OR ' 'E ' FOR ELEC.)" 

42 INPUT T$<3 DISP "TC GET A COMPLETE SET OF DATA. YOU NEE 
D MEASUREMENTS IN BOTH AIR AND WATER. " 

45 DISP "DC AIR FIRST. IN wHAT MEDIUM ARE YOU OPERATING? ( 
ENTER "I" FOR AIR OR '2' FOR WATER)" 

46 INPUT M0 CLEAR 0 BEEP 

47 IF M=2 THEN 49 

48 CS="AIR" 0 GOTO 51 

49 C$= "WATER" 

51 PRINTEF IS 2 0 PLOTTER IS 1 

55 DISP "SET THE DRANETZ FOR ADMITTANCE .SET ON LOWEST SC 
ALE THAT WON'T PEAK DURING RUM.lCCNT)" 

6C PAUSE 

6 5 PRINT IS PRINT 0 PRINT- 
70 PRINT "ADMITTANCE IN ";C$ 

75 GCSUB 7000 

80 IF Mi =2 THEN 50 

85 FG=F6 <s 00=0 0 GOTO ICO 

9 G F 1= F6 0 C1=U 

ICC GCLEAR @ CLEAR 0 BEEP 

105 DISP "SET DRANETZ FOR IMPEDANCE- ( SCAL £ NEEDED ) -HIT C 
CNT " 
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110 PAUSE 

115 PRINT g PRINT g PRINT 
120 PRINT "IMPEDANCE IN " ; C$ 

125 GOSUB 7000 

130 IF M=2 THEN 140 

135 F8=F6 g Q8=Q g GOTO 141 

140 F9=F 6 g (j9=C 

141 CLEAR £ BEEP 

142 PRINTER IS 701 ,76 £ PRINT USING 144 

144 IMAGE 2/ 

145 PRINT " F ( Y -A I R ) = " ; I NT (F 0 ) , , " F ( Y.-WAT) = " ; I NT (F 3 ) , F 
( Z -A I R ) = " ; IN.T ( F8 ) , ," F (Z.-WAT) =" ; INT(F9) 

146 PRINT " C ( Y-AI R ) = " ; CO , w { Y-WATER) = " ;C1 ,,"<*'( Z-AIR) = 
" ;C6 , C(Z- WATER) = ";<jS 

147 PRINTER IS 2 

148 CLEAR @ DISP "THESE ARE THE INITIAL ESTIMATES. (KIT 
CCNT ' WHEN READY TO PROCEED) " 

149 PAUSE 

150 GOSUB 2501 

152 CLEAR £ EEEP 

153 DISP "TO COLLECT ADMITTANCE DATA (FOR ELEC. COUPLING 
) ENTER ' 1' . ENTER "2" FOR IMPEDANCE" 

154 INPUT L2£ PRINTER IS 2 

155 IF L2=2 THEN 158 

156 R$ = "G" £ I $ = " B " £ S5=" ADMITTANCE" (< H$="CONDUCTANCE 
(MI CRCMKOS) " £ G$=" SUSCEPTANCE (MICRGMHCS)" 

157 GOTO 15S 

158 R$="R" £ I$="X" £ S$=" IMPEDANCE" £ G$=" REACTANCE (CH 

MS)" £ H$="RESISTANCE (OHMS)" 

159 IF L2= 1 THEN 176 

16 0 CLEAR @ DISP "SET ON Z AND ENTER SCALE FACTOR TO COL- 
LECT DATA." 

161 INPUT S8£ PRINT S$;" IN ";CS 0 PRINT "SCALE =";£8;" 
OHMS" 

162 CLEAR 

164 GOSUB 8600 

165 IF Y=2 THEN 152 
loo IF M=2 THEN 16 9 

167 R6=£ g X3 = E2 £ XS=E3 £ F4=E4 
16a X6=E 5 £ F5=£6 @ F6=E1 g GOTO 196 

169 R 7=E g X2=£2 § X7 = E3 g F4=E4 

170 X 5=E 5 £ F5=E6 y F9=£l g GOTO 198 

176 CLEAR £ BEEP 

177 DISP "SET ON Y & ENTER SCALE FACTOR IN MICRCMHCS TO 
COLLECT DATA." 

178 INPUT S8e PRINT S?;" IN ";C$ g PRINT "SCALE = " ;S 8;" 
MICRO MHOS" 

179 GOSUB 8600 

ISO IF Y=2 THEN 152 
181 IF M=2 THEN 185 
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Ic2 G6=E Is B 3=E 2 & B8=E3 @ F2=E4 
183 B6=E 5 (s F3=E6 & FG=E 1 0 GCTC 196 
165 G 7=E @ B2 = E2 £ B7=E3 
186 F2=E 4 @ B5 = E5 @ F3=E6 § Fl=El 
156 PRINT USING 199 

199 IMAGE 2/ 

200 BEEP @ DISP "DC YOU DESIFE A LIST OF DA 'I A NEAR RESON 
ANCE? (10G PC INIS) ( L=YES ,2=NC) " 

201 INPUT Y 

202 IF Y=2 THEN 300 
204 IF L2=2 THEN 24G 

210 PRINT "ADMITTANCE DATA FOR " ;C$;" IN MICRO MHOS" £ PR 
INT USING 215 

215 IMAGE 2/ ," FREQUENCY" , 4X ,"REAL" ,6X , "IMAGINARY" 

216 GCTC 265 

24C PRINT "IMPEDANCE DATA FOR " ;C$; " IN OHMS" 0 PRINT US 
ING 215 

265 GCSUB 8100 

300 CLEAR <s BEEP @ DISP "DO YOU DESIRE A PLOT CF DATA? ( 1 
= YES ,2 = NC) " 

305 INPUT Y@ GCLEAR @ CLEAR 
315 IF Y=2 THEN 450 

320 DISP "ENTER: 1= G/R VS. FREw 2= B/X VS. FREG; 3= £ 
VS. G/X VS. R; 4= END PLOTTING LOOP." 

325 INPUT Y 
330 IF Y=2 THEN 365 
335 IF Y=3 THEN 425 
340 IF i = 4 THEN 450 
345 GCSUE 6200 

356 LCIR 0 @ PEN 1 & PENUP 

357 MOVE D ( 1 1 1) + .2* (D (300,1)-D ( 1,1) ) , - ( E* . 8) 

360 LABEL RS;" VS. FREQUENCY IN ";C$ 

361 GOTO 320 
385 GCSUB 8300 

357 LDIR G (9 PEN 1 PENUP 

356 MOVE D(l,l)+.2* (D ( 30u , 1 ) -D ( 1 ,1) ) ,-L-.5*L 

4 00 LABEL 1$;" VS. FREQUENCY IN ";C$ 

401 GCTC 320 

425 GCLEAR & CLEAR S GOSUB 64uG 
431 PEN 1 

442 LDIR C @ PENUP 

443 MOVE - ( . 5*L) ,-L-.2*L 

445 LABEL "INPUT ELECTRICAL " ;SS;" PLOT ( " ; I $ ; " v S . " ; R $ ; " 
) IN ";C$ 

446 GOTO 320 

450 IF L2=2 THEN 645 

451 IF T$="E" THEN 525 

452 GCLEAR g CLEAR 0 DISP "I AM CALCULATING MOTIONAL DAT 
A" 

465 CC=E* .000001/ ( 2*PI * F 6 ) 
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466 LO=l/( ( 2* PI * F6 ) ~ 2*C0 ) 

470 FOR 1= 1 TO 300 

475 D ( I , 2 ) =D ( I ,2).-G 

480 D (I , 3 ) =C (I , 3 ).-D ( I , 1) *B/F6 

485 NEXT I 

4S0 GO SCJB 8700 

500 GOTO 542 

525 GCLEAR 0 CLEAR 0 EISP "I AW CALCULATING MOTIONAL DAT 
A" 

532 C0=8*. 00000 1/(2*PI*F) 

535 FOR 1=1 TO 300 

536 D (I ,2) =C (I , 2).-G 

537 D ( I , 3 ) =D 1 1 , 3 ).-D ( I , 1 ) *3/F 

538 NEXT I 

540 GOSUE 8700 

542 BEEP 0 DISP "DO YOU DESIRE A LIST OF NOTIONAL DATA N 
EAR RESONANCE? ( 1=YES , 2=NO)" 

543 INPUT Y 

545 IF Y=2 THEN 604 

5/5 PRINT "MOTIONAL ADMITTANCE DATA IN ";C$ 0 PRINT USIN 
G 58 0 

580 IMAGE 2/ FREQUENCY" ,4X ."REAL" ,6X /"IMAGINARY" 

555 GCLEAR 0 CLEAR 

6 00 GOSUE 8100 

604 CLEAR 0 BEEP @ EEEP 

6C5 DISF "WAN! A PLOT? ( 1=YES , 2=NC ) " 

610 INPUT Y2 

615 IF Y 2 = 2 THEN 755 

616 PLOTTER IS 705 0 PEN 1 
620 GOSUE 8400 

628 LDIR 0 0 PEN UP 0 PEN 1 

629 MOVE - (. 5*L) ,-L-.2*L 

630 LABEL "MOTIONAL " ;S$;" PLOT FOR ";C$ 

642 GOTO 755 

645 CLEAR 0 GCLEAR 0 DISP "I AM CALCULATING MOTIONAL DAT 

A" 

646 IF T$="M" THEN 656 

647 FOR 1 = 1 TO 300 
646 D( I ,2) =C( I ,2 ) - R 
645 CG=1/(2*PI*F*X) 

650 LG=l/( (2*PI*F6)~2*CG) 

651 D (I , 3) =0 (I ,3).-2*PI*D (I ,1)*LC 

652 NEXT I 

653 GO SUB 8700 

654 C0=AES (CO) 0 L0 =AES ( L 0 ) 

655 EEEP 0 DISP "DC YOU WANT A LIST OF THE DATA NEAR RES 
ON AN CE ? ( 1 = Y ES , 2 =NO ) " 

656 INPUT Y 

657 GOTO 670 

658 L0=X/ ( 2*PI*F6 ) 
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662 FCR 1=1 TO 300 

663 D (I ,2) =D (I , 2).-R 

664 D (I f 3) =D (I , 3).-D (1,1) *2*PI*L0 

6 65 NEXT I 

666 GOSUB 8700 

66 7 £££ P @ DISP " DO YOU WANT A LIST OF DATA NEAR RESCNA 
NCE? ( 1=YES ,2=NQ) ” 

668 INPUT Y 

670 IF Y=2 THEN 706 

68 5 PRINT "MOTIONAL IMPEDANCE DATA IN ";C$ @ PRINT USING 
666 

686 IMAGE 2/ , " FREQUENCY" , 4X , "PEAL " , bX IMAGINARY " 

700 GCLEAR £ CLEAR 

705 GOSUB 8100 

706 CLEAR (§ BEEP @ BEEP 

710 DISP "WANT A PLOT? ( 1 = YES ,2=NC ) " 

715 INPUT Y2 

720 IF Y2=2 THEN 755 

721 PLOTTER IS 70S <§ PEN 1 
725 GOSUB 8400 

737 LDIR 0 i § PENUP 

738 MOVE - (. 5*L) ,-L-.2*L 

740 LABEL "MOTIONAL " ;S$;" PLOT IN " ; C$ 

755 GCLEAR @ CLEAR 

756 PRINT USING 757 

757 IMAGE 3/ 

800 DISP "I AM DOING CALCULATIONS FOR YOU. PLEASE BE PAT 
I ENT." 

SO 3 IF M=2 THEN 1040 
810 IF T $ = " M " THEN 1000 
825 IF L2=2 THEN 660 
640 C 0=E 1/ABS (E4-E6) 
o 45 C1=ABS ( £ 2/ ( 2*PI * F0 ) ) 

6 £ 5 GOTO 88 5 

660 O 8=E 1/AES (£ 4-E6 ) 

665 C1=AES(1/(2*PI*FS*£2) ) 

6 65 K= 1- (F0/F8 ) * 2 
855 Kl=C 3/ (ABS (CO) +C 1) 

S 30 GCSUE 9000 
931 PRINTER IS 2 
940 GOTO 1350 
10U0 K= 1- ( F8/F0 ) "2 
1005 IF L2=2 THEN 1025 
1010 wU=El/ABS (E4-E6) 

1015 GOTO 1 u 30 

1025 i*o = El/AES (E4-E6) 

1030 GOSUB 5200 

1031 PRINTER IS 2 Is GGTG 135u 
1040 IF TS="M" THEN 1150 

lu4 1 BEEP <S DISP "IF YOU HAVE DATA FOR AIR, ENTER 1; ELS 
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E 2" 

1042 INPUT Y 

1043 PAD 

1044 IF L2= 1 THEN 1110 

1045 IF Y=2 THEN 1062 

1060 P6=D4* (D3-D4)/ (C3*R7) 

1061 W 0= F 1 

1062 C I=E 1/ABS ( E 4-E6 ) @ K= 1-(F 1/F9) "2 
1070 GGSUE 9100 

1076 PRINTER IS 2 @ GOTC 1350 
1110 IF Y=2 THEN 1135 
1126 WC=F 1 

1130 P6=D2* (D1-D2 ) / (D1*G7) 

1135 C 1=E 1/AES (E4-E6 ) 

1140 K= 1- (F1/F9) "2 

1145 GGSUB 9100 

1146 PRINTER IS 2 $ GOTO 1350 

1150 RAD 

1151 DISP "IF YCU HAVE DATA FOR AIR, ENTER 1; ELSE 2" 

1152 INPUT Y 

1153 IF L2= 1 THEN 1160 

1154 IF Y=2 THEN 1162 

1155 B=.5*ACS( (R7-R90/D4) 

1156 H6=SIN(B)~2 § H9=COS(B)~2 

1157 P 5= ( ( 1+D3/R0*H 9) ~.5-( L-C3/RC*H8) ~.5)/( ( 1+D3/R0*H9) 

. 5+ ( l-C3/R0*Ha) “ .5) 

1159 P=D4*SIN (2*E) / (4*R9*CS) 0 W0=FS* ( 2*P + ( P '2+ 1) ~ .5 ) 

1160 P6=D4* (D3-D4)/ (C3*R7) 

1 lo2 C5=E1/ABS (E4-E6) £ K= 1- ( F 9/F 1) “2 

1163 PRINT " B= " ; E , " K= " ; K , ” K 1=" ; K 1 , " H8 = " ; rib fi 9= " ; HS 

1164 GCSUB 5300 

1166 PRINTER IS 2 @ GOTC 1350 
1180 RAC 

1182 IF Y = 2 THEN 1210 
1184 E=.5*ACS( (G7-GS)/D2) 

1186 H8=SIN (B) ~ 2 @ H9=CGS(B)~2 

1168 P 5= ( (1+D1/G0*H9) 5-(L-Dl/GC*Ho) 5)/( ( 1+C1/GC*HS) 
. 5+ (1-D1/G0*H&) * .5) 

1190 P=D2*SIN ( 2* E ) / ( 4 *G 9*^ 1 ) <s '.v0= F 1 * ( 2*P+ (P '2 + 1 ) ~ . 5 ) 
1195 P6=C2* (C1-D2)/ (C1*G7) 

1210 w5=El/ABS (E4-E6) @ R = L-(F 9/F 1) "2 

1211 PRINT "£ = ” ;B ,'"K= " ; K t\ 1= " ; K 1 , " H6 = " ;H6 ,"HS=" ;H6 
1215 GCSUB 5300 

1225 PRINTER IS 2 <s GOTO 1350 

1350 DISP "DC YCU DESIRE ANOTHER IN THIS RECluN?(TO GET 
THE OTHER TYPE DATA? ) ( 1=Y ES , 2=i\C ) “ 

1355 INPUT Z1 
1357 PLOTTER IS 1 

1360 IF Z 1=1 THEN 152 

1361 PRINT USING 1362 
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13 62 IMAGE 3/ 

1365 DISP "DC YOCJ DESIRE A RUN IN ANOTHER MEDIUM ?( 1= YES , 
2=NO) " 

1370 INPUT Z1 

1375 IF 21=1 THEN 45 

1360 DISP "JUST FOR THE RECORD, INPUT WATER TEMP, AIR T£ 
MP, TRANSDUCER SER NO, MODEL NO. ," 

1361 DISP "AND TYPE (E OR M ) " 

1385 DISP "INPUT ''O' ' IF INFO IS UNKNOWN" 

13 90 INPUT W,A,S,N,L$ 

1391 PRINTER IS 701,76 

1400 PRINT " aATER TEMP = " ; to , , " AI R TEMP* " ; A , ," SER. NO. = ";S 
, ."MODEL NO.*" ;N , ."TYPE ";L8 
1410 GOTO 9998 

2501 REM *SHUNT* 

2502 CLEAR @ DISP "INPUT A PERCENTAGE OF RESONANT FREQ. 

TO USE TO GET SHUNT (BLOCKED) VALUES." 

2503 DISP "INPUT A VALUE FROM 0 TO 1." 

2504 INPUT El 

2513 CLEAR g BEEF g DISP "WE ARE FINDING SHUNT VALUES FO 
R G/B OR R/X." 

2515 A$="FR" <s B$ = "HZ" 6 S 1= 1 

2519 IF T $ = "M " THEN 2525 

2520 IF M* 2 THEN 2522 

2521 F =F0 *B 1 @ GOTO 25 35 
25 22 F=F 3 *B 1 is GOTO 25 3 5 

2525 IF M* 2 THEN 2527 

2526 F=Fa-61*F6 g GOTO 2535 

2527 F =F 9-B 1*F9 

2535 D$=VAL$ (F.) 

2537 OUTPUT 717 ;A$,D$,D$ 

2536 WAIT 1000 

2540 DISP "SET DRANETZ ERE*. SCALE TO CCvER FREw . ON S YN 
THESIZER .ZERO METERS .SET ON FS.(COMT)" 

2541 PAUSE 

2542 DISP "HERE WE GET THE NORMALIZATION FACTCFS . " 

2543 OUTPUT 7Uy ;"VT1" 

2 548 GCSUB 8549 

2550 G,B,R,X=0 

2551 A $=" FE" £ ES="HZ" t D$=VAL?(F) 

2553 OUTPUT 717 ;A$,C5,E5 

2561 CLEAR Q DISP "SET DRANETZ ON Y ,S£T SCALE, SET ON NCR 
M FOR FILTER + PHASE, PLUG IN TRANSDUCER." 

2562 EEEP (? DISP "ENTER SCALE FACTOR IN MICROPROS . SAmE S 
CALE AS FOR DATA. (^E GET SHUNT G/E)" 

2564 INPUT 3 £ (? CLEAR § DISP "I A w WORKING TO GET SHUNT v 
ALUES" g BEc,P 

2565 ICR 1=1 TO 10 

2566 OUTPUT 717 ;A8,D$,3$ 

2567 r.AIT 1000 
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25bb OUTPUT 709 ?"AC3" 

2569 ENTER 70S ; M(I) 

2570 G=G+N(L) 

2571 OUTPUT 709 ;"AC4" 

2572 ENTER 709 ; N(I) 

2573 B=13+N ( I ) 

2574 NEXT I 

2575 G=G/ ( I- 1 ) t E=E/(I-3) 

2576 G=G*S8/N8 & £=B*SCyN9 @ EEEP 

2576 DISP "SET CRANETZ ON Z; SET SCALE FOR (Aa RESPONSE . k 
OVE TRANS. INPUT. ENTER SCALE FACT OR" 

2579 INPUT Sot CLEAR @ DISP "I kb hOFKING TC GET SHUNT V 
ALUES" 

258G EOF 1=1 TC 10 

2581 OUTPUT 717 ;A$,D$,B$ 

2562 WAIT 1000 

258 3 OUTPUT 709 ;"AC3" 

2 584 ENTER 709 ; K ( I ) 

2585 R=R+M(L) 

2586 OUTPUT 709 ;"AC4" 

2587 ENTER 7u9 ; N(I) 

2588 X=N ( L ) +X 

2589 NEXT I 

2590 R=R/(I-1) @ X = X/(I-3) 

2591 R=R*S8/N8 £ X=X*S8/N9 

2592 IF T $=" M" THEN 2610 

2595 IF to= 1 THEN 2600 

2596 R9=R t X9=X t C9=C t ES-B @ GOTO 26ul 
2600 GC=G t E0=£ t RU=R t XC=X 



2601 


PRINT 


"SHUNT VALUES" 


2602 


PRINT 


USING 2603 


260 3 


I RAGE 


1/ 


2 60 5 


PRINT 


" G0 = " ; GO , , " £ 0 = " ; B 0 , , M R0=";RO, i "XG = ";a0, ,"G9 = " 


•G c 


, " £ 9 = " 


; £9 , ," R9= " ?R 9 , ,"X9 = ";X9 



2606 PRINT t GOTO 2630 

2610 IF S 1=2 THEN 262u 

d. 6 1 1 J 1=G t <j^.= B t J 3 =k t J 4= X 

2612 IF i\=2 THEN 2bl5 

2613 F.=F6+L 1*F 8 § Sl = 2 t GOTO 2550 

2615 F=F9+B 1*F9 t Sl=2 t GOTO 2550 

2620 G= . 5* (G+J 1 ) (s E= . 5* (fi+J 2) t R=.5*(R+U3) & X=.5*(J4+ 

X) t C-CTC 2595 

2630 GCLEAR @ CLEAR t EEEP 

2635 DISP "ENTER''!' IF ALL IS WELL? "2' IF YOU NEED A K 
EPEAT .” 

2o 40 INPUT Y 

2645 IF Y=2 THEN 2501 

2650 RETURN 

7000 REN **F-<*** 

7001 GCLEAR f CLEAR 
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7005 CISP "ENTER LOWER, UPPER FREC. FOR SWEEP. IF RESCRAN 
CE UNKNOWN USE 10000,200000. <FL,F2)" 

7010 INPUT F4,F2 

7011 DISP "ENTER RMS VOLTAGE NEEDED. (30 00 MV FOR LOW FEE 
C. AS LOW AS 2000 MV FOR HIGH)" 

7012 INPUT Afc A 1$=VAL$ (A ) 

7020 CLEAR & BEEP 

7025 A$="AM" 0 B $ = " M R " 

70 35 OUTPUT 717 ;A$,A1$,B$ 

7040 OUTPUT 709 ;"AC3VT3" 

7045 F3=CEIL( (F2-F4)/30C) 

7060 CLEAR @ DISP "I AM WORKING TO GET THE SPECTRUM DAI A 
FOR YOU, .THEN WILL SHOW & MAKE A PIOT" 

7065 FOR 1=1 TO 300 
7070 D ( I , L) =F 4 + 1 * F 3 
7075 A$=" HZ " <e B $ = " F R " 

7083 D$. = VAL$(D (1,1)) 

7085 OUTPUT 717 ;B$,D$,A$ 

7090 OUTPUT 709 ;"AC3VT3" 

7095 WAIT 100 

7100 ENTER 709 ; C(I,2) 

7105 NEXT I 

7110 F5=F4-.l* (F2-F4) @ F6=F 2 + . 1* (F 2-F 4 ) 

7115 F6 = F2-6 . 1* (F2-F4) 

7119 PLOTTER IS 1 g G CL EAR § CLEAR 

7125 SCALE F5 ,F6 ,-. 1 ,1.2 

7130 XAXIS 0 ,2000 ,F4 ,F2 

7135 YAXIS F 4 , .1 ,-.l ,1.2 

7140 FOR 1=1 TO 300 

7145 PEN UP 

7150 PLOT D (I , 1) ,D( I ,2) 

7155 NEXT I 
7158 GRAPH e COPY 

7160 DISP "NEED ANOTHER RuN ? ( 1- YES ,2 = NG ) " 

7161 INPUT Y 

7162 IF 1=1 THEN 7005 

7163 CLEAR Q G CLEAR i? BEEP 

7165 DISP "ENTER DECISION POINT FOR AMPLITUDE (0 TO 1.2) 
. (LESS THAN THE MAX DISPLAYED) ." 

7170 INPUT St CLEAR t BEEP 

7180 PRINT "AMPLITUDE IN vOLTS";" FREQUENCY" 

7 lb 1 PRINT 

7 190 FOR 1=1 IC 30 0 

7195 IF D ( I , 2 ) <S THEN 7210 

7200 PRINT USING 7205 ; D(I ,2) ,D(I ,1) 

7205 IMAGE IX , D . CDDDDCDD , 1 OX .DDDDDC .DC 

7210 NEXT I 

7215 FRINT USING 722U 

7220 IMAGE 3/ 

7225 CLEAR t EEEP 
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7226 DISP "IF K)U ARE READ i TC CONTINUE ENTER "l'; OTHER 
WISE ENTER ' 1 ' . " 

7227 INPUT Y 

7228 IF 1=2 THEN 7160 

7235 DISP "WHAT IS CENTER FREQ, HALFW IQlri TC CONSIDER?" 

7236 INPUT F,N0 CLEAR 0 BEEP 

7240 DISP "I AM WORKING TO GET GOOD FREQ AND Q DATA" 

7241 FOR 1=1 TO 300 

7242 D ( I , 1) =F.-N+I *N/150 
724 3 A$ = "HZ" 9 E $ = " F R " 

7245 D$=VAL$(D (I ,1) ) 

7255 OUTPUT 717 ;E$,D$,A$ 

7260 OUTPUT 709 ;"AC3VT3" 

7265 WAIT 100 

7270 ENTER 70S ; D(I,2) 

7275 NEXT I 

7285 CLEAR 0 G CLEAR 0 DISP "I AM FINDING THE AESOLUTE MA 
X AND FREQ UPPER AND LOWER" 

7 2S0 B 1 ,£4 = C 0 H ,H1 ,H2 ,H3,H4= 1 

7295 B2 , E3= 50 

7300 FOR 1 = 2 TO 300 

7310 IF D (I ,2) <D(H ,2) THEN 7 340 

7315 IF D(I ,2) =D(H ,2) THEN 7325 

7325 A 6 =D ( I , 2 ) 0 F6 = D(I,1) 0 H=I 

7340 NEXT I 

7345 A7=A6./SQR (2) 

7350 FOR 1=1 TO H 
7 355 IF A7 = D ( I ,2) THEN 7405 
7360 IF A7<D ( I ,2) THEN 7385 
7 36 5 IF D ( I , 2) <E ] THEN 7 42 5 
7370 £ 1=C ( I ,2) 0 H 1=1 0 GOTO 7425 
7385 IF E ( I , 2) >E2 THEN 7400 
7 3 SO B2 = D (I ,2) 0 H2=I 
7400 GOTO 7425 

7405 £1 ,£2=D (1,2) 0 Hl,H2=I 9 F7=D(H2,]) 0 GOTO 744U 

7425 NEXT I 

7430 X = (A 7-B 1 ) / (B2-E1) 

74 35 F7=X* ( D ( rl 2 , 1).-C (fil , 1 ) )+D (Hl,l) 

7440 FOR I =H 10 300 
7445 IF A7=D ( I , 2 ) THEN 7495 
7450 IF A7>D( I , 2 ) THEN 747 5 
7455 IF D ( I , 2 ) >B3 THEN 7470 
7460 E 3=D ( I ,2) G H3 = I 
7470 GOTO 7515 

7475 IF D ( I , 2) <=B4 THEN 7490 

7480 £ 4=D (1,2) 0 H4 = I 

7490 GOTO 7515 

7495 B3 ,B4=D (1,2) 9 H3,H4 = I 

7505 L 5=D (H 3 , 1 ) 

7510 GOTO 753C 
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7515 NEXT I 

7520 X= (A7-E4) / (B3-B4) 

7 525 L5=- (X* (D( H4 ,1).-D(H3,1) ))+D (H4,l) 

7530 Q=F6/ ( L5-F7) 

7535 PRINT "CENTER FREQ IS ";F6 

7536 PRINT "Q IS " ;Q 

7538 DISP "IF yCU APE READY TO PROCEED, ENTER ''l'. TO RE 
RUN FOR BETTER VALUES, ENTER " 2 " 

7539 INPUT Y 

7540 IF Y=2 THEN 7235 
7545 RETURN 

8100 REM * * DATA LIST** 

8101 PRINTER IS 2 

8102 FOR 1= 99 TO 199 

6103 PRINT USING 6105 ; D ( I , 1) ,D ( I ,2) ,D ( I , 3 ) 

8105 IMAGE CDCCDD .D ,2X ,D .DDDE , 2X.D.DDCE 

8106 NEXT I 

8107 PRINT USING 8108 
6108 IMAGE 3/ 

6110 RETURN 

6200 REM * * PLOT RE** 

8201 GC LEAR @ CLEAR 

6 202 HS=D ( 1, 1) 15* (D( 3 00 , 1)-D( 1,1) ) 

620 3 H9=D ( 30 0 , 1) + . 1* (D ( 300 , 1 ).-D ( 1,1) ) 

£205 PLOTTER IS 705 @ PEN 1 

6209 SCALE H8 ,H 9 ,- ( . 8 *E ) , E + . 2 5*E 

6210 XAXIS 0 , (D( 300 ,1) -C ( 1, 1) )/lG ,D ( 3 , 3) ,D ( 300 ,1) 

6211 YAXIS D( 1 ,L) ,E/5 ,-(. 6*E) ,E+.2*E 

6212 FEN UP 

8215 FOR 1=1 TO 300 
6216 PLOT D (I ,1) ,D ( I ,2) 

8217 NEXT I 

6216 PENUP £ DEG £ LEIR 0,SIN(90) 

6 221 FOR L1 = D (1,1) TO D(3u0,l) STEP (D ( 30 0 , 1 ) -0 ( 1 , 1) ) / 10 
£ PENUP 

8222 MOVE LI ,- ( . 1S*E) 

822 3 LABEL IN.T(Ll) 

6224 NEXT LI 

6226 LCIR 0 £ PENUP 

8227 MOVE C(l,l) + . 2* ID(300,1)-D (1,1) ) ,-( .7*E) 

6226 LABEL "FREQUENCY (HZ)" 

8230 FENuP £ PEN 1 
d 2 2 2 LDIR 0 

6234 FOR LL=- ( . 8 *E ) TO E+.2*E STEP £/ £ £ PENUP 
62 3 5 MOVE D ( 1 , 1 ) - . 09 * (D ( 3GC , 1 ).-D (1,1)) ,L 1 
6 2 36 LABEL HIT (LI) 

6237 NEXT LI 

6240 LDIR 0,SIN(9Q) 

624 1 MOVE D (1,1).-. 1* (O ( 30u , L)-D ( 1 , 1) ) ,.]*£ 
o24 3 LABEL £8 
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8250 RETURN 
8 200 REN * * P LOT IM** 

8301 GCL EAR § CLEAR 

8 20 3 H6=D( 1 ,1).-. 15* (D ( 300, 1 ) -D ( 1 ,1) ) 

8 204 H9=D ( 300 ,1) + .l*(D(300 f L) -0(1,1)) 

8305 L=MAX (ABS ( E5) ,ABS(E3)) 

8306 PLOTTER IS 705 @ PEN 1 

8312 SCALE H8 ,H9 ,-L- . 5*L ,L+ . 25*L 

8313 XAXIS 0 , (D( 300 ,1)-D (1, 1) )/10 ,D ( 1 , 1) ,D (.200 ,1) 

8214 YAXIS D( 1,1) ,L/5 ,-L-.2*L ,L+.2*L 

8215 PENLP 

8318 FOR 1=1 TO 20 C 
8 225 PLOT D ( I , 1 ) ,D(I,3) 

8320 NEXT I 

8 231 FENUP @ DEG @ LDIR 0,SIN(9G) 

8333 FOR L1=C (1,1) TO 0(200,1) STEP (D ( 20 C , 1 ) -D ( 1 , 1 ) ) / 10 
@ PENUP 

8325 MOVE Ll ,-L- . 2*L 
83 36 LABEL IN.T(L1) 

8337 NEXT Ll 

8340 LDIR 0 (9 PENUP 

8 341 MOVE D (1 , 1) + .2* (D ( 300 , 1)-D ( 1 , 1) ) ,-L-.4*L 
8242 LABEL "FREQUENCY (HZ)" 

8345 PENUP e PEN 1 & LDIR 0 

8 352 FOR Ll=-L- . 2* L TO L-h.2*L STEP L/5 & PENUP 
8 3 55 MOVE D ( 1 , 1) - . 09 * (D (.200 , 1 ) -D ( 1 , 1 ) ) ,L1 
8356 LABEL IN.T(Ll) 

8257 NEXT Ll 

8260 PENUP @ PEN 1 @ LDIR C,SIK(S0) 

8361 MOVE D (1, 1) 1* (D ( 300 , 1 ) -D (1,1) ) ,-(L*.5) 

8268 LABEL G$ 

8380 RETURN 
8 400 REM * *C IRCLES * * 

8401 GCLEAR @ CLEAR 

8402 L 1=MAX ( AES ( E 3) ,AES(E5)) 

8403 L=MAX (ABS (£) ,L1) 

8404 PLOTTER IS 705 0 FEN 1 

8406 SCALE 1. 25* (-L-. 1*L) , 1 . 25* (L+ . 1*L) ,-L- . 2*L,L-h . 15*L 

8407 XAXIS 0 ,L/5 ,-L-. 1*L ,L+. 1*L 
8406 YAXIS 0 ,L/5 ,-L-. 1*L ,L+. 1*L 
8 409 PENUP (5 PEN 1 

8411 FOR 1 = 1 TO 20 0 
8 412 PLOT L (I ,2.) ,D (I , 2) 

8413 NEXT I 

8414 PENUP 

8 4 15 LDIR 0 (§ PEN 1 <s PENUP 
8416 MOVE - (L* . 5) ,-L-. 1*L 
8418 LABEL H$ 

8420 PENUP <3 PEN 1 I? DEG 
8425 LDIR G,SIN(9G) 
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6426 FCR LL=-L-.1*L TO L+.1*L STEP L/ 5 @ PENUP 
8427 MOVE Ll,-L 
6428 LABEL IMT(Ll) 

842S NEXT Ll 

8435 LDIR 0 @ PEN 1 @ PENUP 

8436 FOR LL=-L-.1*L TO L+.1*L STEP L/5 @ PENUP 

8437 MOVE -(.8b6*L),Ll 

8438 LABEL IN.T(L1) 

8439 NEXT Ll 

8445 LDIR G,SIN(9G) 

8446 MOVE - (.95*L) 3*L) 

6448 LAEEL G$ 

6450 RETURN 
8 5 45 REN **NORM* * 

6550 N8 ,N9=0 

6551 FOR 1=1 TO 10 

8557 WAIT 500 

8558 OUTPUT 709 ; "AC 3" 

8559 ENTER 709 ; N(I) 

6560 N8=N8+N ( I) 

6565 OUTPUT 709 ;"AC4” 

6566 ENTER 705 ; M( I) 

6567 NS-N9+MD 

6568 NEXT I 

8569 N8=N8/(I-1) 

8570 N9=-(N9/(I-1) ) 

8572 PRINT USING 8573 
6573 IMAGE 3/ 

8575 RETURN 

8600 REM **vALUES** 

8601 CLEAR @ BEEF 

66 05 IF L2= 1 AND M=2 THEN 86 1C 
8606 IF L2= 1 AND M=1 THEN 6611 

6607 IF L2=2 AND M=1 THEN 8b 12 

6608 IF L2=2 AND to = 2 THEN 6613 

6610 F 6 = F 1 @ & GOTO 8615 

8611 F b=F0 § C=<iC <S GOTO 6615 
86 12 F6 =FS 0 0 GOTO 8615 

6613 F6=F9 0 0=v9 

06 1 5 DISP "PLEASE SELECT THE HALFWIDTH CF INTEREST TO CO 
LLECT DATA. (ENTER HALFWIDTH)” 

6616 INPUT N0 CLEAR 0 DISP "I AM COLLECTING THE DATA." 
6625 FCR 1=1 TC 300 
8629 D (I , 1 ) =F6-N+I *N/ 150 
86 30 A$=" FR" i B $ = " K 2 " 

6632 D$=VAL$(D( I ,1) ) 

6633 OUTPUT 717 ;A$ f D$,E$ 

6634 WAIT 250 

8635 OUTPUT 709 ;"AC3vT3" 

6636 ENTER 70S ; Ho 
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8637 D ( I ,2) =H6*S8/N8 

8638 OUTPUT 709 ;"AC4VT3" 

8639 ENTER 70S ; H9 

8640 D (I ,3) =H9*S8/N9 

8641 NEXT I 

8682 CLEAR g DISP "I AM FINDING KIN/MAX VALUES & ASSCC. 
FREQS'' 

8683 GCSUE 8700 

8684 IF L2=2 THEN 8687 

8685 PRINT "GMAX=";E, FGMX= " ; £ 1 , , " EGKX=" ; E2 , EMX=" ;£ 3 , 
,"FBMX=" ; E 4 , EMI=";E5 , FBMI=" ; E6 



8666 GCTC 8697 

86 67 PRINT "RMAX=" ;E , FRMX=";E1 , ,"XRKX=" ;E2 , XMX = " ;E3 , 



FXMX=" ; E 4 , , " XMI= " ; E 5 , ,"FXMI.= ";E6 

8697 DISP "IF ALL IS *ELL , ENTER ' 1 ENTER '2' TG RETAK 
E DATA." 

8698 INPUT Y 

8699 RETURN 

6 700 REM * *MAXM IN * * 

8703 H ,H 1 , H2= 50 

8705 E=D (50 ,2) @ E 1 , E4 ,E6=D ( 5 0 , 1) @ £ 2 ,E 3 ,E 5 = D 1 5 C , 3) 

8706 FOR 1=51 TC 250 

8710 IF D (I ,2) <E (H,2) THEN 8725 

8711 IF D (I ,2) =D (H ,2) THEN 8720 

8720 £=D(I,2) g E 1 =D ( I , 1 ) & E2=D(I,3) 0 H = I 

8725 IF D (I ,3) <D (Hi ,3) THEN 8740 

8726 IF D (I r 3) =D (HI ,3) THEN 87 35 

6735 E3=D (I , 3 ) 0 E4=D(I,1) g H 1=1 

8740 IF D (I ,3) > = D (H2 , 3) THEN 8751 

3745 H2=I g E5 = D (1,3) 0 E6.=D(I,1) 

8751 NEXT I 

8760 J 1=D (H ,2) 0 J2=(D (H+ 1 , 2) -D ( H-l ,2) )/2 0 J3=(D(H + 1,2) 
+D (H-l ,2).-2*JL)/2 

8752 F6=- ( J 2/ ( 2* J 3 ) ) 

8755 A = J 1+ u 2 +F 6+ J i *F 6 2 

8 770 F6=D (H , 1)+F6* (D(H ,1) -D (H-l ,1) ) 

8775 £ 1=F6 0 E=A 



o 788 D.=ABS ( £ 3-E 5 ) 

6 78 9 IF M= 2 THEN 87 S3 



6 79G IF L2= 1 THEN 8792 



6791 D3=D § GOTO 6796 
3792 D 1 =D 0 GCTC 6796 
o 79 3 IF L2= 1 THEN 8 79 5 



8794 D4= D 0 GCTC 6796 

8795 D2=C 

8 796 PRINT "D 1=" ;Dl , , "D2=";D2 r D3 = " ;D3 , C4 = ";D4 

6799 RETURN 

9000 REM **CALC-A** 

90b 1 PRINTER IS 70 1,76 
90u2 PRINT 0 PRINT 
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9003 PRINT 
S004 PRINT 
SOU PRINT 
;F; "HZ" 

9012 PRINT 

9013 PRINT 

9014 PRINT 

9015 PRINT 

9016 PRINT 



"VALUES MEASURED IN AIR FROM " ;5$;" DATA" 

"ELOCKED CAPACITANCE= ";AES(C0);" FARADS AT " 

"ELECTRICAL QUALITY FACTOR = ";Q0 
"MECHANICAL QUALITY FACTOR = ";<*6 
"ELECTRICAL RESONANCE =";F0;" HZ" 

"MECHANICAL RESONANCE = ";F8;" HZ" 

"DYNAMIC ELECTROMECHANICAL COUPLING COEFFICIE 



NT = 
9 017 
9021 
5031 
9032 
9 050 
9100 



/ 

PRINT "STATIC COUPLING CCEFFICIENT= ";K1 
PRINT " BLOCKED RESISTANCE = " ; R0; "OHMS" 
PRINT USING 9 0 32 
I FACE 2/ 

RETURN 

REM **CALC-W** 



9101 PRINTS 

9102 PRINT 

9103 PRINT 
* DATA" 

9104 PRINT 
Sill PRINT 
" ; F ; " H Z " 
9112 PRINT 
5113 PRINT- 

9114 PRINT 

9115 PRINT 

9116 PRINT 
NT = " ; K 
5119 PRINT 



R IS 701,76 
@ PRINT 

"VALUES FOR MEASUREMENTS IN WATER FROM ";S$; 



"ELOCKED CAPACITANCE = ";A£S(CG);" FARADS AT 

"ELECTRICAL QUALITY FACTOR = ” ; Q ] 

"MECHANICAL QUALITY FACTOR = ";Q$ 

"ELECTRICAL RESONANCE = ";Fl;" HZ" 

"MECHANICAL RESONANCE = " ;F9;" HZ" 

"DYNAMIC ELECTROMECHANICAL COUPLING COEFFICIE 

"FREQUENCY OF OPTIMUM EFFICIENCY = ";WG;" HZ" 



9120 PRINT "MECHANICAL POWER UTILIZATION FACTOR" ; AES ( P 6 ) 



9121 PRINT " STATIC COUPLING CCEFFICIENT-" ; Kl 

9122 PRINT "ELOCKED RESISTANCE =" ; R9 ; "OHMS " 

9132 PRINT USING 9133 

9133 I PAGE 3/ 

9150 RETURN 

S2CU REM **CALC-A/K** 

9201 PRINTER IS 701,76 

9202 PRINT « PRINT 

92C3 PRINT "VALUES MEASURED IN AIR FRCM " ; S 5 ; " DATA" 
9204 PRINT 



9211 PRINT 
;F6 ; "HZ" 

9212 PRINT- 

9213 PRINT 

9214 PRINT 

9215 PRINT 
921b PRINT 



" ELCCKED I NDUCTANCE = ";AES(LC);" HENRIES AT " 

"ELECTRICAL QUALITY FACT CR = " ;c0 
"MECHANICAL QUALITY FACTOR = ";i*6 
"ELECTRICAL RESONANCE =";F0;" HZ" 

"MECHANICAL RESONANCE = ";Fb;" HZ" 

"DYNAMIC ELECTROMECHANICAL COUPLING COLFfICIE 
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NT = " ; K 

5217 PRINT "BLOCKED RESISTANCE® ";R 
9231 PRINT USING S2 32 
5232 INAGE 2/ 

9250 RETURN 

9300 REM * *CALC-W/M* * 

9301 PRINTER IS 701,76 

9302 PRINT 9 PRINT 

9303 PRINT "VALUES FOR MEASUREMENTS IN ft ATE R FROM " ;S$; 
" DATA" 

9304 PRINT 

9311 PRINT "BLOCKED INDUCTANCE = ";A£S(Lu);" HENRIES AT 



";F6; 


"HZ" 






9312 


PRINT 


"ELECTRICAL 


QUALITY FACTOR = " ; w ] 


*313 


PRINT 


"MECHANICAL 


QUALITY FACTOR = " ; 5 


9314 


PRINT- 


" ELECTRICAL 


RESONANCE = " ; F 1 ; " HZ" 


9315 


PRINT 


"MECHANICAL 


RESONANCE = " ; F S ; " H Z " 


9 316 


PRINT 


"DYNAMIC ELECTROMECHANICAL COUPLING - 


NT = 


" ; K 






9318 


PRINT 


"POTENTIAL 


EFFICIENCY = ";F5 


9319 


PRINT 


"FREQUENCY 


CF OPTIMUM EFFICIENCY = " 


9320 


PRINT 


"MECHANICAL 


PCaEF UTILIZATION FACTOR 


9322 


PRINT 


USING 9323 




9 32 3 


IMAGE 


2/ 




9 350 


RETURN 




9998 


DISP " 


.THE END" 




5955 


END 
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